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Abstract 

We study nonlinear dispersive wave systems described by hyperbolic PDE's in M. d 
and difference equations on the lattice The systems involve two small parameters: 
one is the ratio of the slow and the fast time scales, and another one is the ratio of 
the small and the large space scales. We show that a wide class of such systems, 
including nonlinear Schrodinger and Maxwell equations, Fermi-Pasta-Ulam model and 
many other not completely integrable systems, satisfy a superposition principle. The 
principle essentially states that if a nonlinear evolution of a wave starts initially as a sum 
of generic wavepackets (defined as almost monochromatic waves) , then this wave with 
a high accuracy remains a sum of separate wavepacket waves undergoing independent 
nonlinear evolution. The time intervals for which the evolution is considered are long 
enough to observe fully developed nonlinear phenomena for involved wavepackets. In 
particular, our approach provides a simple justification for numerically observed effect 
of almost non-interaction of solitons passing through each other without any recourse 
to the complete integrability. Our analysis does not rely on any ansatz or common 
asymptotic expansions with respect to the two small parameters but it uses rather 
explicit and constructive representation for solutions as functions of the initial data in 
the form of functional analytic series. 

1 Introduction 

The principal object of our studies here is a general nonlinear evolutionary system which 
describes wave propagation in homogeneous media governed either by a hyperbolic PDE's in 
M. d or by a difference equation on the lattice Z d , d = 1,2,3,... is the space dimension. We 
assume the evolution to be governed by the following equation with constant coefficients 

a r U = --L(-iV)U + F(U), U(r,r)| r=0 = h(r), r G M d , (1.1) 

Q 

where (i) U = U(r,r), r G IR d , U G C 2J is a 2J dimensional vector; (ii) L (— iV) is a 
linear self-adjoint differential (pseudodifferential) operator with constant coefficients with 
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the symbol L(k), which is a Hermitian 2J x 2J matrix; (iii) F is a general polynomial 
nonlinearity; (iv) g > is a small parameter. The form of the equation suggests that the 
processes described by it involve two time scales. Since the nonlinearity F (U) is of order 
one, nonlinear effects occur at times r of order one, whereas the natural time scale of linear 
effects, governed by the operator L with the coefficient 1/g, is of order g. Consequently, 
the small parameter g measures the ratio of the slow (nonlinear effects) time scale and the 
fast (linear effects) time scale. A typical example an equation of the form (jl.ljl is nonlinear 
Schrodinger equation (NLS) or a system of NLS. Another one is the Maxwell equation in 
a periodic medium when truncated to a finite number of bands, and more examples are 
discussed below. 

We assume further that the initial data h for the evolution equation (jl.ljl to be the sum 
of a finite number of wavepackets hi, I = 1, . . . , N, i.e. 

h = h 1 + ... + h A r (1.2) 

where the monochromaticity of every wavepacket is characterized by another small pa- 
rameter (3. 

The well known superposition principle is a fundamental property of every linear evolu- 
tionary system, stating that the solution U corresponding to the initial data h as in (J1.2J) 
equals 

U = U x + . . . + U N , for h = h x + . . . + h N , (1.3) 

where U; is the solution to the same linear problem with the initial data h;. 

Evidently the standard superposition principle can not hold exactly as a general principle 
in the presence of a nonlinearity, and, at the first glance, there is no expectation for it to 
hold even approximately. We have discovered though that the superposition principle does 
hold with a high accuracy for general dispersive nonlinear wave systems provided that the 
initial data are a sum of generic wavepackets, and this constitutes the subject of this paper. 
Namely, the superposition principle for nonlinear wave systems states that the solution U 
corresponding to the multi-wavepacket initial data h as in (J1.2J1 equals 

U = Ui + . . . + Uat + D, for h = hi + . . . + h^r, where D is small. 

As to the particular form (jl.ljl we chose to be our primary one, we would like to point out 
that many important classes of problems involving small parameters can be readily reduced 
to the framework of (jl.ljl by a simple rescaling. It can be seen from the following examples. 
First example is a system with a small factor before the nonlinearity 

d t v = -iLv + af (v) , v| t=0 = h, < a < 1, (1.4) 

where initial data are bounded uniformly in a. Such problems are reduced to (jl.ljl by 
the time rescaling r = ta. Note that now g = a and the finite time interval < r < t* 
corresponds to the long time interval < t < r*/ot. 
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Second example is a system with small initial data on a long time interval. The system 
here is given and has no small parameters but the initial data are small, namely 



d t v = — iLv + f (v) , v| t=0 = a h, < a <C 1, where (1.5) 
f (v) = C ] (v) + i1 m+1) (v) + . . . , 

where «o is a small parameter and f*-" 1 ^ (v) is a homogeneous polynomial of degree m > 2. 
After the rescaling v = ctoV we obtain the following equation with a small nonlinearity 



d t V = -iLV + a 



in— l 



ft ] (V) + a f 0(m+1) (V) + . . . , V| t=0 = h, (1.6) 



which is of the form of (jl.4|) with a = a m 1 . Introducing the slow time variable r = tot™ 1 
we get from the above an equation of the form (jl.lj) . namely 

W = "^T LV + [ f (m) ( v ) + «of (m+1) (V) + . . .] , V| t=0 = h, (1.7) 



where the nonlinearity does not vanish as cto — > 0. In this case g = at™ 1 and the finite time 
interval < r < corresponds to the long time interval < t < l^-i with small «o < 1. 

Very often in theoretical studies of equations of the form (jl.ljl or ones reducible to it a 
functional dependence between g and (3 is imposed, resulting in a single small parameter. 
The most common scaling is g = (3 2 . The nonlinear evolution of wavepackets for a variety 
of equations which can be reduced to the form (jl.ljl was studied in numerous physical and 
mathematical papers, mostly by asymptotic expansions of solutions with respect to a single 
small parameter similar to 0, see [HI, d, PI, 120], E3, EB], E3, El, EH, Ei> H3 
and references therein. Often the asymptotic expansions are based on a specific ansatz 
prescribing a certain form to the solution. In our studies here we do not use asymptotic 
expansions with respect to a small parameter and do not prescribe a specific form to the 
solution, but we impose conditions on the initial data requiring it to be a wavepacket or a 
linear combination of wavepackets. Since we want to establish a general property of a wide 
class of systems, we apply a general enough dynamical approach. There is a number of general 
approaches developed for the studies of high- dimensional and infinite-dimensional nonlinear 
evolutionary systems of hyperbolic type, [TO], [TJ], E3, E3> EU> E3> E3> EH> S3) 
j45j ) and references therein. We develop here an approach which allows to exploit specific 
properties of a certain class of initial data, namely wavepackets and their linear combinaions, 
which comply with the symmetries of equations. Such a class of the initial data is obviously 
lesser than all possible initial data. One of the key mathematical tools developed here for 
the nonlinear studies is a refined implicit function theorem (Theorem I4.25J) . This theorem 
provides a constructive and rather explicit representation of the solution to an abstract 
nonlinear equation in a Banach space as a certain functional series. The representation is 
explicit enough to prove the superposition principle and is general enough to carry out the 
studies of the problem without imposing restrictions on dimension of the problem, structural 
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restrictions on nonlinearities or a functional dependence between the two small parameters 
g,(3. 

As we have already stated the superposition principle holds with high accuracy for linear 
combinations of wavepackets. A wavepacket h (/3, r) can be most easily described in terms 
of its Fourier transform h (/3, k). Simply speaking, wavepacket h (/3, k) is a function which is 
localized in /3-neighborhood of a given wavevector k* (the wavepacket center) and as a vector 
is an eigenfunction of the matrix L(k), details of the definition of the wavepacket can be 
found in the following Section 2. The simplest example of a wavepacket is a function of the 
form 

h 09, k) = £T d h (^J^) Sn (k.) , k e M. d , (1.8) 

where g n (k*) is an eigenvector of the matrix L (k*) and h (k) is a Schwartz function (i.e. it is 
infinitely smooth and rapidly decaying one). Note that the inverse Fourier transform h (/3, r) 
of h (j3, k) has the form 

h (p, v) = h (/5r) e ik * r g„ (k.) ,reR d , (1.9) 

where h (r) is a Schwartz function, and obviously has a large spatial extension of order 
We study the nonlinear evolution equation (jl.lj) on a finite time interval 

< r < r*, where r* > is a fixed number (1-10) 

which may depend on the L°° norm of the initial data h but, importantly, r* does not depend 
on g. We consider classes of initial data such that wave evolution governed by is 
significantly nonlinear on time interval [0, r*] and the effect of the nonlinearity F (U) does 
not vanish as g — > 0. We assume that j3, g satisfy 

2 

< (3 < 1, < q < 1, — < Ci with some Ci > 0. (1.11) 

e 

The above condition on the dispersion parameter £- ensures that the dispersive effects are 
not dominant and do not suppress nonlinear effects, see [7 j for a discussion. 

To formulate the superposition principle more precisely we introduce first the solution 
operator S (h) (r) : h — ► U (r) which relates to the initial data h of the nonlinear evolution 
equation (jl.lj) the solution U (t) of this equation. Suppose that the initial state is a multi- 
wavepacket, namely h = h^, with h;, / = 1, . . . , N being "generic" wavepackets. Then for 
all times < r < r* the following superposition principle holds 

||D (r)|| B = sup ||D (r)|| LOO < C 5 ^tj for any small 5 > 0. (1.13) 

0<r<r* j3 

Obviously, the right-hand side of (j!.13j) may be small only if g < C\f3. There are examples 
(see 0) in which D (r) is not small for g = C\f3. In what follows we refer to a linear 
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combination of wavepackets as a multi-wavepacket, and to wavepackets which constitutes the 
multi-wavepacket as component wavepackets. 

The superposition principle implies, in particular, that in the process of nonlinear evo- 
lution every single wavepacket propagates almost independently of other wavepackets even 
though they may "collide" in physical space for a certain period of time and the exact solution 
equals the sum of particular single wavepacket solutions with a high precision. In particu- 
lar, the dynamics of a solution with multi-wavepacket initial data is reduced to dynamics of 
separate solutions with single wavepacket data. Note that the nonlinear evolution of a single 
wavepacket solution for many problems is studied in detail, namely it is well approximated 
by its own nonlinear Schrodinger equation (NLS), see [H], j23|, EH], GUI, EH|, ED], EH, 
and references therein. 

The superposition principle (|1.12|) . (|1.13|) can also be looked at as a form of separation of 
variables. Such a form of separation of variables is different from usual complete integrability, 
and its important factor is the continuity of spectrum of the linear component of the system. 
The approximate superposition principle imposes certain restrictions on dynamics which 
differ from usual constraints imposed by the conserved quantities as in completely integrable 
systems as well as from topological constraints related to invariant tori as in KAM theory. 

Now we present an elementary physical argument justifying the superposition principle. 
If nonlinearity is absent, the superposition principle holds exactly and any deviation from it is 
due to the nonlinear interactions between wavepackets, so we need to estimate their impact. 
Suppose that initially at time r = the spatial extension s of every composite wavepacket is 
characterized by the parameter (3~ l as in (jl.9j) . Assume also (and it is quite an assumption) 
that the component wavepackets during the nonlinear evolution maintain somehow their 
wavepacket identity, group velocities and spatial extension. Then, consequently, the spatial 
extension of every component wavepacket is propositional to and its group velocity v j is 
proportional to p -1 . The difference Av between any two different component group velocities 
is also proportional to g -1 . The time when two different component wavepackets overlap in 
space is proportional to s/ \Av | and, hence, to q/(3. Since the nonlinear term is of order 
one, the magnitude of the impact of the nonlinearity during this time interval should be 
proportional to g/0, which results in the same order of magnitude of D. This conclusion is 
in agreement with the estimate of magnitude of D in (11.13)) (if we set 5 = 0). 

The rigorous proof of the superposition principle we present in this paper is not based 
on the above argument since it implicitly relies on a superposition principle in the form 
of an assumption that component wavepackets can somehow maintain their identity, group 
velocities and spatial extension during nonlinear evolution which by no means is obvious. In 
fact, the question if a wavepacket or a multi-wavepacket structure can be preserved during 
nonlinear evolution is important and interesting question on its own right. The answer 
to it under natural conditions is affirmative as we have shown in [7j. Namely, if initially 
solution was a multi-wavepacket at r = 0, it remains a multi-wavepacket at r > 0, and every 
component wavepacket maintains its identity. Therefore a wavepacket can be interpreted as 
a quasi-particle which maintains its identity and can interact with other quasi-particles. This 
property holds also in the situation when there are stronger nonlinear interactions between 
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wavepacket components which do not allow the superposition principle to hold, see [7] for 
details. 

The proof we present here is based on general algebraic-functional considerations. The 
strategy of our proof is as follows. First, we prove that the operator S (h) in (jl.l2|) is 
analytical, i.e. it can be written in the form of a convergent series 

S (h) = V °° S U) (h j ) , h 3 = h, . . . , h (j copies of h), 

where (hP) is a j-linear operator applied to h. Now we substitute h in with the sum 
of hi as in (jl.2j) . Considering for simplicity the case N = 2 and using the polylinearity of 
<S^ we get 

((h a + h 2 ) 2 ) = ((hO 2 ) + 2S® (h!h a ) + S (2) ((h 2 ) 2 ) , . . . , 
implying after the summation 

S (h) = ((h^ 2 ) + S (3) ((h x ) 3 ) + . . . + S (2) ((h 2 ) 2 ) + ((h 2 ) 3 ) + . . . s cr 
= «S(hi)+,S(h 2 )+«S cr , 

where S cr is a sum of all cross-terms such as (hih 2 ) etc. The main part of the proof is 
to show that every term in S cr is small. An important step for that is based on the refined 
implicit function theorem ( Theorem I4.25|) which allows to represent the operators in the 
form of a sum of certain composition monomials, which, in turn, have a relatively simple 
oscillatory integral representation. Importantly, the relevant oscillatory integrals involve the 
known initial data h rather than unknown solution U. The analysis of the oscillatory in- 
tegrals shows that there are two mechanisms responsible for the smallness of the integrals. 
The first one is time averaging, and the second one is based on large group velocities (in 
the slow time scale) of wavepackets. Remarkably, if wavepackets satisfy proper genericity 
conditions, every cross term is small due one of the above mentioned two mechanisms. Im- 
portantly, the both mechanism are instrumental for the smallness of terms in S cr , and the 
time averaging alone is not sufficient. We obtain estimates on terms in S cr which ultimately 
yield the estimate (jl.l3|) . Since the smallness of interactions between waves under nonlinear 
evolution stems from high frequency oscillations in time and space of functions involved in 
the interaction integrals, we can interpret it as a result of the destructive wave interference. 
The above sketch shows that the mathematical tools we use in our studies are (i) the theory 
of analytic functions and corresponding series of infinite-dimensional (Banach) variable, and 
(ii) the theory of oscillatory integrals. 

We would like to point out that the estimate (|1.13|) for the remainder in the superposition 
principle is quite accurate. For example, when the estimate is applied to the sine- Gordon 
equation with bimodal initial data, it yields essentially optimal estimates for the magnitude 
of the interaction of counterpropagating waves. These estimates are more accurate than 
ones obtained by the well known ansatz method as in [SH], and the comparative analysis is 
provided below in Example 1, Section 2.2. 

To summarize the above analysis we list important ingredients of our approach. 
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• The spectrum of the underlying linear problem is continuous. 

• The wave nonlinear evolution is analyzed based on the modal decomposition with re- 
spect to the linear component of the system because there is no exchange of energy 
between modes by linear mechanisms. Wavepacket definition is based on the modal 
expansion determining, in particular, its the spatial extension and the group velocity. 

• The problem involves two small parameters (3 and g respectively in the initial data 
and coefficients of the equations. These parameters scale respectively (i) the range of 
wavevectors involved in its modal composition, with (3~ scaling its spatial extension, 
and (ii) g scaling the ratio of the slow and the fast time scales. We make no assumption 
on the functional dependence between /3 and g, which are essentially independent and 
are subject only to inequalities. 

• The nonlinear evolution is studied for a finite time r* which may depend on, say, 
the amplitude of the initial excitation, and, importantly, is long enough to observe 
appreciable nonlinear phenomena which are not vanishingly small. The superposition 
principle can be extended to longer time intervals up to blow-up time or even infinity if 
relevant uniform in f3 and g estimates of solutions in appropriate norms are available. 

• Two fast wave processes (in the chosen slow time scale) attributed to the linear operator 
L and having typical time scale of order g can be identified as responsible for the 
essential independence of wavepackets: (i) fast time oscillations which lead to time 
averaging; (ii) fast wavepacket propagation with large group velocities produce effective 
weakening of interactions which are not subjected to time averaging. 

The rest of the paper is organized as follows. In the following Section 2 we formulate ex- 
act conditions and theorems for lattice equations and partial differential equations and give 
examples. In Section 3 we recast the original evolution equation in a convenient reduced form 
allowing, in particular, to construct a representation of the solution in a form of convergent 
functional operator series explicitly involving the equation nonlinear term. In Section 4 we 
provide the detailed analysis of function- analytic series used to get a constructive represen- 
tation of the solution. Section 5 is devoted to the analysis of certain oscillatory integrals 
which are terms of the series representing the solution. Note that when making estimations 
we use the same letter C for different constants in different statements. Finally, the proofs 
of Theorems 12 . 151 and 12 . 1 91 are provided in Section 6. more examples and generalizations are 
given in Section 7. For reader's convenience we provide a list of notations in the end of the 
paper. 

2 Statement of results 

In this section we consider two classes of problems: lattice equations and partial differential 
equations. After Fourier transform they can be written in the modal form which is essentially 



7 



the same in both cases. We formulate the exact conditions on the modal equations and present 
the main theorems on the superposition principle. We also give examples of equations to 
which the general theorems apply, in particular Fermi-Pasta-Ulam system and Nonlinear 
Schrodinger equation. 

2.1 Main definitions, statements and examples for the lattice equa- 
tion 

The first class of evolutionary systems we consider involves systems of equations describing 
coupled nonlinear oscillators on a lattice Z d , namely the following lattice system of ordinary 
differential equations (ODE's) with respect to time 

<9 T U (m, r) = --LU (m, r) + F (U) (m, r) , U (m, 0) = h (m) , m G Z d , (2.1) 
£? 

where L is a linear operator, F is a nonlinear operator and g > is a small parameter (see 
jH]). To analyze the evolution equation (|2.1|) it is instrumental to recast it in the modal 
form (the wavevector domain), in other words, to apply to it the lattice Fourier transform 
as defined by the formula 

U (k) = ^ U ( m ) e ~ im ' k > where k e I' 71 ' A d > ( 2 - 2 ) 

mGZ d 

k is called a wave vector. We assume that the Fourier transformation of the original lattice 
evolutionary equation (|2.1j) is of the form 

<9 T U (k, r) = --L (k) U (k, r) + F (tj) (k, r) ; U (k, 0) = h (k) for r = 0. (2.3) 

Here, U (k, r) is 2 J- component vector, L (k) is a k-dependent 2 J x 2 J matrix that cor- 
responds to the linear operator L and F ^U^ is a nonlinear operator, which we describe 

later. The matrix L (k) and the coefficients of the nonlinear operator F ^U^ in (J2.3)) are 2n- 

periodic functions of k and for that reason we assume that k belongs to the torus M. d / (2n'L) d 
which we denote by [—tt, n] d . The k-dependent matrix L (k) determines the linear operator 
L and plays an important role in the analysis. We refer to L (k) as to the linear symbol. 
Since ()2.3|) describes evolution of the Fourier modes of the solution, we call (j2.3j) modal 
evolution equation. 

We study the modal evolution equation (|2.3|) on a finite time interval 

< r < r„ (2.4) 

where r* > is a fixed number which, as we will see, may depend on the magnitude of the 
initial data. The time t* does not depend on small parameters, it is of order one and is 
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determined by norms of operators and initial data; it is almost optimal for general F since 
there are examples when r* is of the same order as the blow up time of solutions. To make 
formulas and estimates simpler we assume without loss of generality that 

r* < 1. (2.5) 

For a number of reasons the modal form (J2.3)) of the evolution equation is much more suitable 
for nonlinear analysis than the original evolution equation (J2.1J) . This is why from now on 
we consider the modal form of evolution equation h2.ty for the modal components U (k, r) as 
our primary evolution equation. 

First, as an illustration, let us look at the simplest nontrivial example of (|2.3|) with J = 1 
corresponding to two-component vector fields on the lattice Z d . A two-component vector 
function U (m) of a discrete argument m G Z d has the form 



U(m) 



17+ (m) 
U_(m) 



m G Z d . (2.6) 



In this example L (k) in ()2.3|) is a 2 x 2 matrix, and we assume that for almost all k it has 
two different real eigenvalues oj- (k) and uj + (k) (the dependence of u± (k) on k is called the 
dispersion relation) satisfying the relation uo- (k) = —uo + (k), namely, 

L (k) g c (k) = u ( (k) g ( (k) , u c (k) = (u (k) , C = ±, (2.7) 

where, evidently, g^ (k) are the eigenvectors of L (k). These eigenvalues (k), ( = ±, are 
27r-periodic real valued functions 

uj ( (kt + 27T, k 2 , . . . , k d ) = . . . = u c (ki, k 2 ,...,k d + 2tt) = u c (k t , k 2 , . . . , k d ) . (2.8) 

The simplest nonlinearity in ()2.3|) is a quadratic nonlinear operator F ^U^ = F^ (jj 2 ^j 
which is given by the following convolution integral 



F^ (uA) (k) = J ( k , k) (Ux (kO U 2 (k")) dk', 



(2.9) 



k'e[-7r,7r] d ; k'+k" 



where k = (k, k"), ^k, kj is a quadratic tensor (susceptibility) which acts on vectors 

Ui,U2- We refer to the case J = 1 as the one-band case since the corresponding linear 
operator is described by a single function uo (k). 

A particular example of ()2.3|) is obtained as a Fourier transform of the following Fermi- 
Pasta- Ulam equation (FPU) (see |T2], jSZ], jSl) describing a nonlinear system of coupled 
oscillators: 

d T x n = - (y n - y n _i) , (2.10) 
Q 

d T y n = ~ {Xn+1 - »n) + "2 (X n +1 ~ ^nf + " 3 (x n +l ~ %nf , Tl G Z. 

Q 
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Note that an equivalent form of (j2.1U|) (with a 2 = 0) is the second order equation 



1 



g 2 q 



X , 



{x n X<n— l) ) • 



(2.11) 



In this example d — 1, k = k and elementary computations show that the Fourier transform 
of the FPU equation (|2.10|) has the form of the modal evolution equation ([2.3)1 . ()2.9|) where 



U 



x 

y 



iL (k) 







(l - e" ifc ) 



(l-e" ifc )* 



u c (k) = 2C 



k 



X 



(2) (k, k', k") Ui {k') U 2 (k" 



a 2 1 



-ik J 



1 



-ijfc' 



sm ■ 







(2.12) 



5i (A;') x 2 (fc") 



and a similar formula for x^ 3 -* (see (17.5)) ). 

Now let us consider the general multi-component vector case with J > 1 which we refer 
to as J -band case for which the system ()2.3|) has 2J components, and instead of ()2.7|) we 
assume that L (k) has eigenvalues and eigenvectors as follows: 



L (k) g n , c (k) = u n>c (k) g njC (k) , u niC (k) = C^ n (k) , C = ±, n—l,...,J, 



(2.13) 



where u n (k) are real- valued, continuous for all k functions, and eigenvectors g n ^ (k) 6 C 2J 
have unit length in the standard Euclidean norm. We also suppose that the eigenvalues are 
numbered so that 

uJn+i (k) > co n (k) > 0, n = 1, . . . , J - 1, (2.14) 

and we call n the 6and index. Note that the presence of ( = ± reflects a symmetry of the 
system allowing it, in particular, to have real- valued solutions. Such a symmetry of dispersion 
relation u n (k) occurs in photonic crystals and many other physical problems. 
Note that ([2.13)1 implies that the following symmetry relation hold: 



u n _ c (k) = -u nX (k) , n = 1, . . . , J. 
We also always assume that the following inversion symmetry holds: 
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-k) 



(k). 



(2.15) 



(2.16) 



Remark 2.1 Assuming \2. Inl and \2. 1 6|) we suppose that the dispersion relations uo^ (k) 
have the same symmetry properties as the dispersion relations of Maxwell equations in periodic 
media, see ^-J^, J3J^. We would like to stress that these symmetry conditions are not imposed 
for technical reasons but because they are consequences of fundamental symmetries of physical 
media. Such symmetries arise in many problems including, for instance, the Fermi-Pasta- 
Ulam equation, or when L (k) originates from a Hamiltonian H (p, q) = \ (Hi (p 2 ) ) + ^H 2 (q 2 ). 
In the opposite case if it is assumed that \2. lty and 112.1b}) never hold, the results of this paper 
hold and the proofs, in fact, are simpler. The case with the symmetry is more difficult and 
delicate because of a possibility of resonant nonlinear interactions. 
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There are values of k for which inequalities (J2.14j) turn into equalities, these points require 
special treatment. 

Definition 2.2 (band-crossing points) We call k a band-crossing point if u n+ i (k ) = 
uj n (k ) for some n or U\ (k ) = and denote the set of band-crossing points by a. 

Everywhere in this paper we assume that the following condition is satisfied. 

Condition 2.3 The set o~ of band-crossing points is a closed nowhere dense set in M. d with 
zero Lebesgue measure, the entries of the matrixL (k) are infinitely smooth functions o/k ^ a 
and uj n (k) are continuous functions o/k for all k and are infinitely smooth when k ^ a. 

Observe that for k ^ o all the eigenvalues of the matrix L (k) are different and the 
corresponding eigenvectors g n ^ (k) of L (k) can be locally defined as smooth functions of 
k ^ a as long as L (k) is smooth. 

Remark 2.4 The band- crossing points are discussed in more details in £7]/, Jlf. Here we only 
note that generically the singular set a is a manifold of the dimension d — 2, see l^). A 
simple example of a band-crossing point is k = in \2.12]) . 

Since we do not assume the matrix L (k) to be Hermitian, we impose the following 
condition on its eigenf unctions which guarantees its uniform diagonalization. 

Condition 2.5 We assume that the 2J x 2J matrix formed by the eigenvectors g n ^ (k) of 
L (k), namely, 

E (k) = [g li+ (k) , gl ,_ (k) , . . . , g Ji+ (k) , g Ji+ (k)] 
is uniformly bounded together with its inverse 

sup ||S(k)|| , sup 1 1 ST 1 (k)|| < C= for some constant C~,. (2-17) 

Here and everywhere we use the standard Euclidean norm in C 2J . 

Note that if the matrix L (k) is Hermitian for every k, the eigenvectors form an orthonor- 
mal system. Then the matrix H, which diagonalizes L, is unitary and (J2.17j) is satisfied with 
Cg = 1. Everywhere throughout the paper we assume that Condition \2.5\ is satisfied. 

We introduce for vectors u e C 2J their expansion with respect to the basis g n ^: 

j j 

* ( k ) = E E <c w 8n.c ( k ) = E E ( k ) > ( 2 - 18 ) 

n=l C=± n=l C=± 

and we refer to it as the modal decomposition of u(k), and call the coefficients u n ^ (k) the 
modal coefficients of u(k). In this expansion we assign to every n, ( a linear projection 
U n £ (k) in C 2J corresponding to (k), namely 

n n>c (k) u (k) = u niC (k) g„ i( (k) = u n , c (k) , n = 1, . . . , J, C = ±- (2-19) 
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Note that these projections may be not orthogonal if L (k) is not Hermitian. Evidently the 
projections n n ( j (k) are determined by the matrix L (k) and therefore do not depend on the 
choice of the basis g n ^ (k). Projections U n ^ (k) depend smoothly on k a (note that we 
do not assume that the basis elements g n ^ (k) are defined globally as smooth functions for 
all k ^ cr, in fact band-crossing points may be branching points for eigenf unctions, see for 
example Ij.) They are also uniformly bounded thanks to Condition 12.51 



1/2 



<C 3 |V|, V G C 2J , k£o 



[2.20) 



We would like to point out that most of the quantities are defined outside of the singular set 
a of band-crossing points. It is sufficient since we consider U (k) as an element of the space 
Li of Lebesgue integrable functions and the set a has zero Lebesgue measure. 

The class of nonlinearities F in (J2.3j) which we consider can be described as follows. F is 
a general polynomial nonlinearity of the form 



mp 

p (p) = p {m) (^ m ) ' with mp - 2 ' 



(2.21] 



m=2 



where m-linear operators F^ are represented by integral convolution formulas similar to 
(J22J), namely 

(U 1; . . . , U m ) (k, r)= f X (m) (k, k) d (k') . . . U m (k^ (k, £)) d^" 1 ^, (2.22) 

3 m = [-7r,4 m - 1)d , (2.23) 
~S m ~ l ) d k = dk' ... dk (m ~ 1} (2.24) 

(2.25) 



where the domain 
and we use notation 



and 



(27i) (m - 1)d 

k {m) (k, kj = k - k' - . . . - k (m - 1} , k = (k', . . . , k (m) ) 



Condition 2.6 (nonlinearity regularity) The nonlinear operator F ( U ) defined bu \2.21\) 
satisfy 



X 



(m) I 



(2tt 



■^33 sup \\x {m) {k,k',...M m) )\\<C x , ro = 2,3, 



(2.26) 



k,k' ...,k( m ) 



where, without loss of generality, we can assume that C x > 1. The norm x^ (k, kj of the 
tensor x^ with a fixed k as a m-linear operator from (C 2J ) m into (C 2J ) is defined by 

(2.27) 



X {m) (k, k) 


= sup 


X {m) (m)(X!,. 


• ; ^-m) 




l*jl<l 
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where as always, \ ■ \ stands for the standard Euclidean norm. The tensors x^ (k, kj are 

assumed to be smooth functions o/k, k', . . . , k*" 1 ) ^ a, namely for every compact K C W 1 \ o 
and for all m = 2,3, ... . 

\V l X {m) (k, k', . . . , k< m >) | < Ckjl if k, k', . . . , k< m > 6 K, I = 1, 2, . . . , (2.28) 

where t/ie vector composed of all partial derivatives of order I of all components of 

the tensor x with respect to the variables k, k', . . . , k.( m \ 

/,From now on all the nonlinear operators we consider are assumed to satisfy the nonlin- 
earity regularity Condition 12.61 

Remark 2.7 At first sight, since g is a small parameter, one might think that the linear term 
in \2. 1)) with the factor - is dominant. But it is not that simple. Indeed, since all eigenvalues 

o/L(k) are purely imaginary the magnitude of e~~e L ^h (k) which represents the solution 
of a linear equation (with F = 0) is bounded uniformly in g. A nonlinearity F alters the 
solution for a bounded time which is not small for small g. Therefore the influence of 
the nonlinearity can be significant. This phenomenon can be illustrated by the following toy 
model. Let us consider the partial differential equation for a scalar function y (x, r) : 

d T y = —d x y + y 2 , y {x, 0) = h (x) . 
Q 

Its solution is of the form 

y(x,r) = V / ' v (2-29) 

1-rhlx-Z) 

and regularly it exists only for a finite time. The solution \2. 2ty) shows that the large coef- 
ficient - enters it so that the corresponding wave moves faster with the velocity - along the 
x-axis but the wave's shape does not depend on g at all. For the NLS with the initial data 
h (k) = h (k, (3), g = 1 , and the coefficient - at the linear part, the nonlinearity balances 
the effect of dispersion leading to emergence of solitons, see ^ for a discussion. 

To formulate our results we introduce a Banach space E = C ([0, r*] ,Li) of functions 
v (k, r), < r < r*, with the norm 

||v(k,r)|| E =||v(k,r)|| c([art]iLi) = sup f |v(k,r)|dk. (2.30) 

0<T<T, J[-7T,7r] <i 

Here Li is the Lebesgue function space with the standard norm defined by the formula 

|v(k)|dk. (2.31) 

id 
7T,7Tj 

The following theorem guarantees the existence and the uniqueness of a solution to the modal 
evolution equation ()2.3)1 on a time interval which does not depend on g (see Theorem 15 . 41 for 
details). 




13 



Theorem 2.8 (existence and uniqueness) Let the modal evolution equation \2.'J\) satisfy 
the Condition \2.5l and let h G L lt h < R. Then there exists a unique solution U = 

Li 

Q (lij of H~3J) which belongs to C 1 ([0,r*] , L%). The number r* > depends on R, C x and 
Cs and it does not depend on g. 

Now we would like to formulate the main result of this paper, a theorem on the super- 
position principle, showing that the generic wavepackets evolve almost independently for the 
case of lattice equations. To do that, first, we define an important concept of wavepacket. 

Definition 2.9 (wavepacket) A function h (/3, k) which depends on a parameter < (3 < 
1, is called a wavepacket with a center k* if it satisfies the following conditions: 



(i) It is bounded in L\ uniformly in (3, i.e. 

h(/3, 



<C h , 



Li 



(2.32) 



(ii) It is composed of modes from essentially a single band n, namely for any < e < 1 
there is a constant C e > such that 

h (k) - h_ (k) - h + (k) < ca h c (k) = n n , c h (k) , c = ±, (2.33) 

Li 

and h^ (/3, k) is essentially supported in a small vicinity of^k*, where k* is the wavepacket 
center, namely 



Ik-ck.l^ 1 



h c (0, k) 



dk < C e (3. 



(2.34) 



(Hi) The wavepacket center 'k* is not a band-crossing point, that is k*/ ^ a, and the following 
regularity condition holds: 



V k h c (0, k) 



dk < C e p 



-l-e 



'Ik-Ck.l^ 1 -^ 

In the above conditions (ii) and (Hi) C e does not depend on [3, < /3 < 1. 



(2.35) 



The simplest example of a wavepacket in the sense of Definition 12.91 is a function of the 
form 



h ( ((3, k) = (3- d h 



k-Ck, 



,C (k) , C = ±, 



(2.36) 



where (k) is a Schwartz function, that is an infinitely smooth, rapidly decaying function. 
Another typical and natural example of a wavepacket h centered at k* is readily provided by 



h (/3, k) = n n , + (k) h ,+ 03, k) + n n ,_ (k) h 0i _ k) 



(2.37) 
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where h ^ (/3, k) is the lattice Fourier transform of the following function 

h 0iC (m, (3) = e* k * m <I> c ((3m - r ) g, C = ±, (2.38) 

where g is a vector in C 2J , projection Tl n £ is as in (|2.19|) with some n, vectors m, ro G M. d 
and (r) being an arbitrary Schwartz function (see Lemma f7.2|) . 

Our special interest is in the waves that are finite sums of wavepackets and we refer to 
them as multi-wavepackets. 

Definition 2.10 (multi-wavepacket) A function h(f3, k), < (3 < 1, is called a multi- 
wavepacket if it is a finite sum of wavepackets h; as defined in Definition \2.tA namely 

h(/3,k) =£}h, ((3,k), (2.39) 
i=i 

and we call the set {k*i} of all the centers k*/ of involved wavepackets center set ofh. 

In what follows we will be interested in generic multi-wavepackets such that their centers 
are generic. The exact meaning of this is provided below in the following conditions. 

Condition 2.11 (non-zero frequency) We assume that every center k^ of a wavepacket 
satisfies the following condition 

u ni (K,)^Q, l = l,... 7 N h . (2.40) 

Condition 2.12 (group velocity) We assume that all centers k^ ; / = l,...,Nh, of the 
multi-wavepacket h as defined in Definition \2.1(A are not band-crossing points, and the gra- 
dients VkW ni . (k*^) (called group velocities) at these points satisfy the following condition 

Vk^n^ (k*zj - VkU ni2 (k*/ 2 ) ^ when h ^ l 2 , (2.41) 
indicating that the group velocities are different. 

We also want the functions (dispersion relations) uj ni (k) to be non-degenerate in the sense 
that they are not exactly linear, below we give exact conditions. 
Consider the following equation for n and 9 

6u ni {K)-Cu> n {8K) = 0, C = ±l, (2.42) 
where the admissible 6 have the form 



Til 



9 = ^C (i) , C (i) = ±1, m < m F , (2.43) 
i=i 

mp is the same as in (|2.21|) . In the case when in the series (j2.21j) some terms vanish, 
we take in (J2.43|) only m corresponding to non-zero F^ m \ 
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Condition 2.13 (non-degeneracy) Given a point k* = and band n\ we assume that 
dispersion relations u n (k) are such that all solutions n,9 of \2.^2^ are necessarily of the form 

n = n h 9 = (. (2.44) 

Definition 2.14 (generic multi-wavepackets) A multi-wavepacket h as defined in Def- 
inition Y2.1U\ is called generic if the centers k^, I = l,...,Nh, of all wavepackets satisfy 
Conditions \2.11\ and \2. 12\ and the dispersion relations u n (k) at every k*; and band ni satisfy 
Condition HH51 

We introduce now the solution operator Q mapping the initial data h into the solution 

< R 

owing 



U = Q (iij of the modal evolution equation (|2.3|) ; this operator is defined for 
according to Theorem 12.81 The main result of this paper for the lattice case is the fo 



statement. 



Theorem 2.15 (superposition principle for lattice equations) Suppose that the ini- 
tial data h of A2.3\) is a multi-wavepacket of the form 

N h 

h= Vh ; , A^max hi < R, (2.45) 



1=1 



satisfying Definition \2.1(A where h is generic in the sense of Definition \2. l\ Let us assume 
that 

B 2 1 1 

— < C, with some C, 0<(3<-, 0<g<-. (2.46) 

Then the solution U = Q (h^ to the evolution equation \2. S\) satisfies the following approxi- 
mate superposition principle 

\Z=1 / 1=1 

with a small remainder D (r) satisfying the following estimate 

Dfr 



sup 

0<T<T* 



<C * |ln0|, (2.48) 



where e is the same as in Definition \2.iA and can be arbitrary small, r* does not depend on 
3, g and e. 

The most common case when ()2.46|) holds is g = ft 2 , a discussion of different scalings is 
provided in and [7j. 

Observe that solutions to the original evolution equation ()2.1|) with the initial data ()2.39|) . 
()2.H8|) satisfy the superposition principle if the wave vectors k*; in (j2.H8j) satisfy ()2.41j) . ()2.42|1 
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and are Schwartz functions. It turns out, that the evolution of every coefficient u n ^ (k) 
of the solution as defined by (J2.18|) can be accurately approximated by a solution a relevant 
Nonlinear Schrodinger equation (NLS), see [2H]- Therefore Theorem l2. 15l provides a reduction 
of multi-wavepacket problem to several single-wavepacket problems. 

We also would like to stress that though (3 is small the nonlinear effects are not small. 
Namely, there can be a significant difference between solutions of a nonlinear and the corre- 
sponding linear (with F (U) being set zero) equations with the same initial data for times 

T = T*. 

Recall that up to now we analyzed the nonlinear evolution in the modal form ()2.3|) for 
U (k, r). To make a statement on the nonlinear evolution for the original evolution equation 
1)2. ip . i.e. in terms of the quantities U (m, r), we introduce U (h) (m) as the inverse Fourier 

transform of the solution Q (h^j (k) of the modal evolution equation (J2.3)) . Recall that the 

inverse Fourier transform corresponding to ()2.2|) is given by the formula 

U (m) = (2n)~ d I d e imk U (k) dk, (2.49) 

J [— 7T,7r] 

and when applying the inverse Fourier transform we get back the original lattice system ()2.ip 
from its modal form (|2.3p . The convolution form of the nonlinearity makes the lattice system 
invariant with respect to translations on the lattice Z d . Using Theorem 12.151 and applying 
the inverse Fourier transform together with the inequality 

||U|| Loo < (2ny d U (2.50) 

we obtain the following statement. 

Corollary 2.16 Let the evolution equation \2. 1)) be obtained as the lattice Fourier transform 
of \2.3\l . If h is given by \2. Sty) where every $^ (r) is a Schwartz function (that is an 
infinitely smooth, rapidly decaying function) then U (h) is a solution to the evolution equation 
\2.1)) . Ifh = hi + . . -+h.N h and every hi is given by \2. Sty) then the approximate superposition 
principle holds: 

U(h) =U(h 1 ) + ... + U(h iV J+D, (2.51) 
with a small coupling remainder D (r) satisfying 

sup ||D(r)|| Loo <^-^, (2.52) 

0<T<T„ P 

where 5 > can be taken arbitrary small. 

As an application of Theorem 12. 151 let us consider the Fermi- Pasta-Ulam equation ()2.1()j) . 
We impose the initial condition for ()2.1()p 



rih n h 

x n (0) = ^oi {fin - n) e^ n + cc, y n (0) = £ ifin - n) e ik * ;ri + cc, n G Z, (2.53) 

i=i 
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where \I/oz (r) , (r) are arbitrary Schwartz functions, and r; are arbitrary real numbers, cc 
means complex conjugate to the preceding terms and assume that g, (3 satisfy (|2.46|) . For 
any given k*i there are two eigenvectors g± (k*i) of the matrix L (A;*/) in ()2.12|) given by ()7.3|) 
and corresponding terms in ()2.53|) can be written as 



^01 



[$-,ig- ( k *i) + $+,/g+ (k*i))e 



In this case all requirements of Definition 12.101 are fulfilled, and ()2.53j) defines a multi- 
wavepacket. Note that the multiwavepacket ()2.53|) involves Nh = 2nh wavepackets with 
2nh wavepacket centers tik^i,® = ±. To satisfy Condition 12. 121 the wavepacket centers k*i 
must satisfy 



cos ■ 



sin 



cos ■ 



sm 



if/^j. 



(2.54) 



To check if the centers k*i satisfy Condition 12.131 we consider the equation 

3 



sm 



<<*/ 



-c 



sm z 



k*i 



0, z 



i=i 



c 



0) 



±i. 



(2.55) 



Evidently the possible values of z are —3,-1,1,3. Since the equation 3|sin</>| = |sin(30)| 
has the only solution <fi = on [0,7r/2] the equation ()2.55|) has the only solution z = (. 
Consequently, all points k*i ^ satisfy Condition I2.13| and Theorem 12.151 applies. The 
initial data for a single wavepacket solution have the form 



x&,n,i (0) 

y#,n,i (o) 



(ftn - n) {hi) + cc, n e z, ■& = ±. 



(2.56) 



According to this theorem and Corollary 12.161 the solution to ()2.10j) . ()2.53|) equals the sum 
of solutions of (|2.1()j) with single wavepacket initial data, that is 



x n (r) = ^ x #> n > 1 ^ + Dl ' n ( r ) ' Vn ( r ) = S V{> ^ 1 ( r ) + D ^ n ^ 

0=± 1=1 i?=± 1=1 

where D n is a small remainder satisfying 

Q 



sup sup [\D ljU (r) 

0<t<t* n 



D 2 , n (r) \<C 6 



1+6 



(2.57) 



(2.58 



with arbitrarily small positive S. Hence, the following statement holds. 



Theorem 2.17 (superposition for Fermi-Pasta-Ulam equation) If every $>i^ (r) is a 
Schwartz function, and the wavevectors k*i ^ satisfy \2. 54 ), then the solution x n (r) , y n (r) 
of the initial value problem for the Fermi-Pasta- Ulam equation \2.10)) with multi-wavepacket 
initial condition \2.5§) is a linear superposition of solutions x n> i (r) , y n)i (r) of the same equa- 
tion with single-wavepacket initial condition 112.56]) up to a small coupling term D\^ n (r) , £>2,n ( r ) 
satisfying fi?. 57}) , h2. 5fyl with arbitrary small 5 > and r* which do not depend on (3, g,S. 
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Note that solutions x$^i (r) with different I resemble 2n^ solitons which originate at 
different points r ; and propagate with different group velocities. According to ()2.57j) . (|2.58|) 
all these soliton-like wavepackets pass through one another with very little interaction, see 
Fig. Note that Theorem 12 . 1 51 shows that this phenomenon is robust in the class of general 




Figure 1: In this picture two wavepackets are shown with different "centers" and k^- 
The values of fc*i and k* 2 are proportional to the frequences of spatial oscillations. Though 
the wavepackets overlap in physical space, they pass one through another in the process of 
nonlinear evolution almost without interaction if their group velocities are different. 

difference equations on the lattice Z, and that it persists under polynomial perturbations of 
the nonlinearity as well as perturbations of the linear part of the equation (|2.11j) as long as 
they leave the linear difference operator non-positive and self-adjoint. Observe also that the 
evolution of every single wavepacket is nonlinear, and it is well-approximated by a properly 
constructed NLS (we intend to write a proof of this statement for general lattice systems in 
another article; see (231 f° r a particular case). For example, for a special choice of tyji the 
solution x nj i (r) can be well approximated by a soliton solution of a corresponding NLS. 

2.2 Main statements and examples for semilinear systems of hy- 
perbolic PDE 

In this subsection we consider nonlinear evolution equation involving partial differential (and 
pseudodifferential) operators with respect to spatial variables with constant coefficients in 
the entire space M. d . There is a great deal of similarity between such nonlinear evolution 
PDE and the lattice nonlinear evolution equations considered in the previous section. In 
particular, we study first not the original PDE but its Fourier transform, modal evolution 
equation, and the results concerning the original PDE are obtained by applying the inverse 
Fourier transform. 

Recall that for functions U (r) from L\ (M d ) the Fourier transform and its inverse are 
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defined by the formulas 

U (k) = ! U (r) e~ irk dr, where k G R d , (2.59) 



U (r) = / U (k) e irk dr, where r G R d . (2.60) 

(27f) J]Rd 

Similarly to ()2.3j) we introduce the following modal evolution equation 

d T U (k, r) = -~L (k) U (k, t) + F (ft) (k, r) , U (k, 0) = h (k) , k G R d , (2.61) 

where (i) U (k, r) is a 2 J-component vector-function of k, r, (ii) L (k) is a 2 J x 2 J matrix 

function of k, and (iii) F ^Ijj is the nonlinearity. We assume that the 2 J x 2 J matrix L (k), 

k G R d , has exactly 2 J eigenvectors g n ^ (k) with corresponding 2 J real eigenvalues u) n £ (k) 
satisfying the relation (j2~T3]l . flZHP , (f2T5jl . (f2~Tfift . (pHTjl . We also assume the matrix L (k), 
k G M d , to satisfy the polynomial bound 

|L(k)| <C(l + |k| p ). (2.62) 

The singular set a for L (k) is as in Definition 12.31 with the only difference that functions 
u n £ (k) are defined over R d rather than the torus [— 7r,7r] d , and, consequently they are not 

periodic. The nonlinearity F (U) has a form entirely similar to (j2.21|) : 



^ (p) = ^ (m) (^ m ) ' ( 2 - 63 ) 

m=2 

with F^ being m-linear operators with the following representation similar to ()2.22|) : 

(u a , . . . , U m ) (k) = I X {m) (k, fc) Uj (k') . . . U m (k( m ) (k, fc)) d^" 1 ^, (2.64) 

t/ ID) m 



where k.^ ys.,kj is defined by the convolution equation (|2.25|) . d is defined by ()2.24|) and 
D m in (ETMI) is now defined not by (I2~23|) but by 

D m = M (m ~ 1)rf . (2.65) 

The difference with (12. 3|) now is that the involved functions of k, k' etc., are not 27r-periodic, 
D m in ()2.64|) is defined by ()2.65|) instead of ()2.23|) . and the tensors fk, k) satisfy the non- 
linear regularity Condition 12.61 without the periodicity assumption. The functions (k^J 
in (12.64)) are assumed to be from the space L\ = Li (R d ^j with the norm 



U 



|v(k)|dk. (2.66) 
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We seek solutions to (|2.61|) in the space C 1 ([0, r*] , L\) with < r* < 1. 

Applying the inverse Fourier transform to the modal evolution equation (|2.61|) we obtain 
a hyperbolic 2J-component systems in W d of the form 



T V (r, r) = L (-iV r ) U (r, r) + F (U) (r, r) , U (r, 0) = h (r) 

Q 



(2.67) 



Note that since L (k) satisfies the polynomial bound 12.62)) we can define the action of the 
operator L (— iV r ) on any Schwartz function Y (r) by the formula 



-iV r )Y(k) = L(k)Y(k) 



(2.68) 



where, in view of (j2.62J) . the order of L does not exceed p. If all the entries of L (k) are poly- 
nomials, such a definition coincides with the common definition of the action of a differential 
operator L(— iV r ). In this case L (— iV r ) defined by ()2.68|) is a differential operator with 
constant coefficients of order not greater than p. 

The properties of the modal evolution equation (|2.61|) are completely similar to its lattice 
counterpart and are as follows. The existence and uniqueness theorem is similar to Theorem 



Theorem 2.18 (existence and uniqueness) Let equation \2. 61)) satisfy conditions ( 2.1''/\j 



< R. 



and (OOP and h e L x = L x (R d ), 

modal evolution equation 112.61)) in the functional space C 1 
pends on R, C x and C-=. 



Then there exists a unique solution to the 
, tJjLx). The number r* de- 



Here is the main result for the semilinear hyperbolic systems of PDE which is completely 
similar to Theorem 12.151 

Theorem 2.19 (principle of superposition for PDE systems) Let the initial data of 
the modal evolution equation \2. 61)) be a multi-wavepacket, i.e. the sum of Nh wavepackets 
hi as in \2.4ty satisfying Definitions \2. yi \2.1(A Suppose that g, (3 satisfy condition jj2.46\) . 

Assume also thath is generic in the sense of Definition \2. 1J\ Then the solution U = Q ^hj to 

the modal evolution equation \2. 61)) satisfies the approximate linear superposition principle, 
namely 

/Nh \ Nh 



^2.69) 



. i=i 



with a small remainder D (t) 



sup 

0<T<T„ 



DM 



i=i 



< C t 

Li ~ p 



Q 



l+e 



Ml 



(2.70) 
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where e is the same as in Definition \2. 91 r* does not depend on (3,g and e. The solutions 
U (h) (r, r) of the space evolution equation \2. 6'7| ) are obtained as the inverse Fourier trans- 
form of Q ( h ) and they satisfy the approximate linear superposition principle, namely 



U (h) = U (hi) + . . . + U (h;vj + D, (2.71) 
with a small coupling remainder D (r) satisfying 

snp \\D{T)\\ Loo <C e -^- e \\nP\, (2.72) 

0<T<r» p 

where e > is the same as in Definition \2.9i and can be arbitrary small. 

Example 1: Sine-Gordon and Klein-Gordon equations with small initial data Let 

us consider the sine-Gordon equation (see j2H|) 

8 2 t u = d 2 u - sinw (2.73) 

with small initial data 

u (r, 0) = Pb , d t u (r, 0) = [3h, < 1. (2.74) 

First, we recast this the equation into our framework by rescaling the variables 

u = pu u (3 2 t = t. (2.75) 

Since sm.f3U\ = f3U\ — ^(3 3 Uf + P 5 f(Ui), where evidently f (Ui) is an enitire function, we 
can recast the equation (J2.73|) into the following form 

%Ui = j, [d 2 x U x -U,]+j 2 [qUf + (3 2 f (170] • (2-76) 

1 /2 

We introduce then a linear pseudodifferential operator A = (I — d 2 ) with the symbol 

1/2 j^^rm. 

(1 + k 2 ) and rewrite the equation ()2.76|) as the following system 

d T Ut = ^AU 2 , d T U 2 = -i^i + A- 1 [qUl + (3 2 f (t/x)] , (2.77) 
with the initial data 

I7i(0) =ho, U 2 (0) = h 1 , (2.78) 
where ho and hi are assumed to be of the form 

z (r, 0) = h , p (r, 0) = hi, h 3 = £ *ii (P r ~ r i) + cc ' 3 = 0, 1, (2.79) 

i=i 
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in one-dimentional case with r = r, k = k. Evidently, the relations with the initial data of 
(gZZSD are 

b = h , h = Ah x . 
Notice that the system (|2.77|) is of the form (|2.b7|) with 



q = /3 2 , LU = 
F (U) = A- 1 



AU 2 
-AU X 



qUf 



F(U)=F (U)+/? 2 F 1 (U) 
Ft (U) = A- 1 



(2.80) 





f{Ui) 



Observe now that L has only one spectral band with the dispersion relation and eigenvectors 
given by 



-vd 
1 



7? = ±1, 



^01 



u (k) = (I + k 2 f\ &, (k) = = 2- 1 ' 2 
and there is no band-crossing points. We use expansion in the basis g± 

[$ + , zg+ + $_, z g_]e ik - r . (2.81) 

to represent initial data ()2.78|) and (|2.79|) . The equation ()2.42|) takes here the form 

(l + ^) 1/2 A = C(l + A 2 ^) 1/2 ; C = ±L 

Obviously, this equation has only solutions A = ( and Condition 12.131 is fulfilled. Condition 
I2~T21 holds if 

^^7rS> for '^' ort >^' (2 - 82) 



which is equivalent to 



7^ k*i for I' ^ I, and k*i ^ for all I. 



^2.83) 



Equation (|2.77|) can be written in the integral form (|3.3|) with mp = oo and by Theorem 
15.41 it has unique solution U for r < r*. If we replace F (U) in (|2.8U|) by F (U) we obtain 



d T Vy = ^AV 2 , d T V 2 = -j 2 AV x + A-'qV^ 

where we take the initial data to be as in ()2.78|) . namely 

V 1 (0)^h o , Vt{0) = hi. 



(2.84) 



(2.85) 



Equations (|2.84|1 can be obtained by replacing sinw in ()2.73|) by the cubic polynomial u — u 3 /6 
producing the quasilinear Klein-Gordon equation (see [HE]). Observe that the solutions to 
the sine-Gordon and the Klein-Gordon equations with small initial data are very close. To 
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see that, note that the operator / (U) (k) is bounded in L\ for U (k) which are bounded in 



L\. Therefore the norm of the neglected term is small, namely 



P 2 f (U) 



by Remark 14. 8| the solutions of (|2.77jl and (|2.84jl are close, namely 

ll^i-^illi +\\U 2 -V 2 \\ L <C(3 2 , 0<r<r, 

II II Ijfv) II II Uoo ' 



< C{3 2 . Thus, 



(2.86) 



According to Theorem 12. 191 the superposition principle is applicable to the equation ()2.84|) 
with initial data as in (|2.85|) . and the following statements hold. 

Theorem 2.20 (Superposition for Klein-Gordon) Assume that the initial data ho, hi 
in \2.85)) are as in \2. 71]) . Then the solution {Vi, V 2 } to the system \2.84\) satisfies the linear 
superposition principle, namely 



V x (r, r) = E V W ( r ' T ) + Dl ( r > r ) ' V * ( r ' T ) = E E V W ( r ' r ) + ° 2 ( r > T ) ' ( 2 - 87 ) 



0=± z=i 



=± i=i 



where {Vj.,^ (r, r) , (r, r)} is a solution to \2. 84}) with the one-wavepacket initial condi- 
tion 

Vi^i (r, 0) (R s ik r 

^2,^ (r, UJ 

where $^ (r) are arbitrary Schwartz functions. If \2.8ty holds, the coupling terms Di,D 2 
satisfy the bound 



(2.* 



sup [HDiCr) 



|D 2 (r)|| L J<C^ = ^- 5 , 



(2.89) 



0<T<T* 

where r* and C' s do not depend on (3, and 5 can be taken arbitrary small. 

Using ()2.86|) we obtain a similar superposition theorem for the sine-Gordon equation. 

Theorem 2.21 (Superposition for sine-Gordon) Assume that the initial data ho, hi in 
(PTTgp are as in $F7$) . Then the solution {Ui,U 2 } to PTTTp , JE75| ) satisfies the linear 
superposition principle, namely 



£/! (r, r) = EE r ) + Dl ( r ' r ) ' ^ ( r ' r ) = E E ^.i ( r > r ) + °2 (r, r) , 

0=± Z=l j?=± 2=1 

where £/i,#,z (r, r) , f/ 2 .i?,z ( r ? r ) a solution of \2. 77| ) i/ie one-wavepacket initial condition 



Ui A i (r, 0) 



(/3r-r0g^ r + cc,^ = ± 



where (r) are arbitrary Schwartz functions. If A2.8ty) holds, the coupling terms D 1; D 2 
satisfy the bound i2.8!A) . 
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Note that a theorem completely similar to Theorem I2.2UI holds also for a generalized 
Klein-Gordon equation where qVf is replaced by an arbitrary polynomial P (Vi). Hence, the 
superposition principle holds for the sine-Gordon equation ()2.73|) with a small initial data and 
a strongly perturbed nonlinearity as, for example, when sinu is replaced by sinu + /5 _1 u 4 + 

/rV. 

We would like to compare now our results and methods with that of [38] where the 
interaction of counterpropagating waves is studied by the ansatz method. Pierce and Wayne 
considered in [3H] the sine-Gordon equation in the case of small initial data which have the 
form of a bimodal wavepacket. In our notation it corresponds to the case when g = (3 2 , 
rih = 1 in ()2.79|) . when two wavepackets, corresponding to ■& = + and ■& = — , have exactly 
opposite group velocities. They proved that the bimodal wavepacket data generate two waves 
which are described by two uncoupled nonlinear Schrodinger equations with a small error. 
The magnitude of the error given in jSH] (which we formulate here for the solution U\ of 
the rescaled equation ()2.76|) ) is estimated by C/3 1 ^ 2 on the time interval < r < r ( 
or < t < tq/3^ 2 ). Note that our general Theorem 12. 191 when applied to the special case 
of the sine-Gordon equation ()2.76|) provides a better estimate of the coupling error, namely 
Cg/f3 1+s = C(3 l ~ 5 in fl£Hj| with arbitrary small 5, for the same time interval. Notice that 
the estimate (J2.72|) given in Theorem 12.191 is almost optimal, since it is possible to construct 
examples when the coupling error is greater than c(3 1+s with arbitrary small S. 

We would like to point out that the general mechanism responsible for the wavepacket 
decoupling is the destructive wave interference, this mechanism is subtle though general. 
We treat the destructive wave interference by taking into account explicitly all nonlinear 
interactions of high-frequency waves. In our approach we use the exact representation of a 
general solution in the form of a functional-analytic operator monomial series, every term 
of the series is explicitly given as a multilinear oscillatory integral operator applied to the 
initial data. A key advantage of such an approach is that it allows to estimate wavepacket 
coupling as a sum of contributions of highly oscillatory terms and to get a precise estimate of 
magnitude of every term. In contrast, the well known "ansatz" approach as, for instance, in 
[3S] and [22], requires to find a clever ansatz with consequent estimations of the "residuum" 
in an appropriate norm. Our approach can naturally treat general tensorial polynomial 
nonlinearities F of arbitrary large degree Np and any number of wavepackets, whereas finding 
a good ansatz which allows to estimate the residuum in such a general situation would be 
difficult. For readers interested in detailed features of one- wavepacket solutions to the sine- 
Gordon equations, we refer to [32] , [HE] and [55] . 

Example 2: Nonlinear Schrodinger equation. The Nonlinear Schrodinger equation 
(NLS) with d spatial variables ([12], [IE], [IS]) has the form 

d T z (r, r) = i-7 (-iV) z (r, r) + a \z\ 2 z (r, r) , z (r, 0) = h (r) , re R d , (2.90) 
Q 
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where a is a complex constant, 7 (— iV) is a second-order differential operator, its symbol 
7 (k) is a real, symmetric quadratic form 

7 (k) = 7 (k, k) = lijkkj, 7 (~iV) z = - Y lijdrAjZ- 

To put the NLS into the framework of this paper we introduce the following two-component 
system 

d T z + (r, r) = i-7 (-iV) z + (r, r) + az.z 2 + (r, r) , (2.91) 
Q 

d T z_ (r, r) = -i-7 (iV) z_ (r, r) + a*z + z 2 _ (r, r) , 
Q 

z + (r, 0) = h (r) , z_ (r, 0) = h* (r) , r G R d , 

where a* denotes complex conjugate to a. Obviously if z (r, r) is a solution of (|2.9U|) then 
z + (r, r) = z (r, r), z_ (r, r) = (r, r) gives a solution of (|2.91|) . Using the Fourier transform 
we get from (|2.90|) 



d T z (k, r) = i-7 (k) 5 (k, r) + a(z*z2) (k, r) , k G 
' 



^2.92) 



Now the band-crossine set a = {k G M d : 7 (k) = 0}. We assume that the quadratic form 
7 is not identically zero. The Fourier transform of (|2.91|) takes the form of (J2.67|) with 



U 



U- 



L(k)U 

(k) = | 7 (k)|,F( 3 ) (U 



7 (k) 




a[z + (U 

a* (z_ (U ) z_ ( U ) z . ( U 










7 (-k) _ 




_ U- _ 


)4E 





To satisfy the requirements of Condition 12. 141 we have to take the wave vectors k*/ ^ a so 
that 

2 7(k*0 , , 2 7 (k,,/) , 
7 (k*j, •) 7^ — — — \~j~7 (k*i', •) if I t 1 1 



l7(k*z)| l7(k*i')| 
which provides (|2.41j) . Since 

| 7 A - C |7 (Akrt)l = |7 (Ml [A - C |A| 2 ] , 

and A is odd, every point k^ ^ a satisfies Condition 12.131 If the quadratic form 7 is not 
singular, that is det7 7^ 0, then condition (|2.93p . which ensures that group velocities of 
wavepackets are different, holds when 

7 ( k *0 . / 7 (k*jQ , 

It Ml l7( k *r)l 

In this case Theorem 12.191 is applicable, and generic wavepacket solutions of the NLS are 
linearly superposed and propagate almost independently with coupling O {(3). More precisely, 
as a corollary of Theorem 12. 191 we obtain the following statement. 
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Theorem 2.22 (Superposition for NLS) Assume that initial data of the NLS \2.9U\) have 
the form h = h\ + . . . + h,N h 

h (r) = e ik - m $ z , + (/3r - r ) + e - ik - m $^ (/3r - r ) , I = 1, . . . , N h 

where (r) are arbitrary Schwartz functions. Assume also that det 7 ^ and the vectors 
k*/ satisfy conditions 

TW^O, l = l,...,N h ; K^K, if i^v. 
Then solution z = z (h) is a linear superposition 

z(h) = z (Tii) + ... + Z (h Nh ) + D 

with a small coupling term D 

sup \\D{r)\\ L ( wd) <Cs-^, 

0<t<t» v ' p 

where 5 > can be taken arbitrary small. 



We note in conclusion, that the superposition principle reduces dynamics of multi-wavepacket 
solutions to dynamics of single-wavepacket solutions; we do not study dynamics of single- 
wavepacket solutions in this paper. Note that the theory of NLS-type approximations of 
one-wavepacket solutions of hyperbolic PDE is well-developed, see [21], jSOI, |H] jUIl, !S], 
and references therein. Relevance of different group velocities of wavepackets for smallness 
of their interaction was noted in |29j . 



2.3 Generalizations 

Note that in a degenerate case when the function u ni (k) is linear in the direction of k* the 
equation ()2.42|) for C, = 1 has many solutions for which 6 7^ ±1 and Condition 12. 131 does not 
hold. It turns out, that if Condition 12. 13l for dispersion relations u n (k) at k* is not satisfied, 
still we can prove our results under the following alternative condition. We consider here the 
case of PDE in the entire space M. d and k e M. d . 

Condition 2.23 (complete degeneracy) The series \2.21\) has only with odd m. 

The wavevectors k^ and functions u) ni (k), 1 = 1,..., N h , have the following three properties: 

(i) There exists 5 > such that for every l\ 7^ I2, the following inequality holds: 

V k u; n!i {viKh) - V k w„ i2 (z/ 2 k*« 2 ) > 5, (2.94) 

for any odd integers v\,v 2 = 1,3,.... 

(ii) There exists 5 > such that z/k*; does not get in a 5 -neighborhood of a for any odd 
integer v and any 1 = 1,..., Nh. 
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(Hi) For any positive integer odd number 9 and any k^, for any n the following identities 
hold: 

V k uj n (6Ki) = V k w„ (Ki) , (2.95) 



u n (0k,,) = 6u n (Ki) . (2.96) 

A nontrivial examples, where the above Condition 12.231 is satisfied, is given below. 
We give here a generalization of Definition 12.141 

Definition 2.24 (generic multi-wavepackets) A multi-wavepacket h as defined in Defi- 
nition WU\ is called generic if (i) the centers k^ ; 1 = 1,..., Nh, of all wavepackets satisfy 
Conditions \2.11\ and \2.12i - (ii) either the dispersion relations u n (k) at every k*, and band n\ 
satisfy Condition \2.1?A or they satisfy Condition \2.23l 

The statement of Theorem l2. 19l remains true if Condition l2.14l is replaced by less restrictive 
Condition 12.241 namely the following theorem holds. 

Theorem 2.25 Let the initial data of the modal evolution equation \2. 61)) be a multi-wavepacket, 
i.e. the sum of wavepackets h; as in \2.4^j) satisfying Definitions \2. iA [2.1 (A Suppose that 
\2.46\l holds. Assume also that h is generic in the sense of Definition \2. 2~$\ Then the so- 
lution U = Q [h^ to the modal evolution equation \2. 61]) satisfies the approximate linear 
superposition principle, namely \2.6iA) . \2. 7(Jj) , \2. 71]) and \2. 72]) hold. 



The proofs we give in this paper directly apply to more general Theorem 12.251 
Another generalization concerns the possibility to shift independently initial wavepackets. 
If initial data involve parameters as in ()2.79j) it is possible to prove that C e in (|2.48J) . ()2.7()j) 
and 1)2.72)1 does not depend on G M. d if the functions are Schwartz functions. Most 
of the proofs remain the same, but several statements have to be modified, and we present 
proofs in a subsequent paper. 

One more generalization concerns the smoothness of initial data. It is possible to take 
initial data h/ (r) with a finite smoothness rather than from Schwartz class. Namely, consider 
weighted spaces L 1(l with the norm 

||v|| Lia = / (l + |k|r|v(k)| dk, a>0. (2.97) 

Obviously, large a corresponds to high smoothness of the inverse Fourier transform v (r) . 
Then if functions (k) have the form ()2.36)) with h\ (k) = hi^ (k) from the class L\ A the 
inequality ()2.70)) can be replaced by 



sup 

0<T<T* 



D(r) <C £ -^|ln/3|+a/? s , (2.98) 

Li jj 



28 



where s > and e > have to satisfy restriction - < a. This generalization requires minor 
modifications in the proofs and in conditions (J2.33)) and ()2.34j) C e /3 has to be replaced by 
C t (3 s . In particular, if a = 1, g = (3 2 and s — 1/2 the right-hand side of ()2.98|) can be 
estimated by C^/? 1 ' 2-61 with arbitrary small e±. 

More generalizations which involve the structure of equations are discussed in Sections 
7.3 and 7.4. Now we give an example where Condition 12.231 is applicable. 



Example 3: Semilinear wave equation. Let us consider a semilinear wave equation 
with d spatial variables 



d 2 T z (r, r) = - Az (r, r) + -d Xl z 3 (r, r) , r G R d , 



g* 



Q 



(2.99) 



where A is the Laplace operator, a is an arbitrary complex constant, g = (3 . We introduce 
the operator A = y/—A which is defined in terms of the Fourier transform, it has symbol |k|. 
We rewrite ()2.99|) in the form of a first-order system 



d T z (r, r) = -Ap (r, r) , r G R d ; 


d T p (r, r) = — Az (r, r) + aA^d^z 3 (r, r) . 
Q 



(2.100) 



The linear operator A x d Xl has the symbol -j^p h is a zero order operator. We rewrite 
()2. 100ft in the form of piTTj) where 



U 



z 
P 



, -iL(-iV r )U 






A ' 




z 








-A 









•'( 




y = a 








. P . 





-A-id Xl z 3 



Using the Fourier transform we get (|2.fjlj) with 



U 



z 
V 



, -iL(k)U 
1 



Ikl 



z 
V 



#(3) fu 



-iafci 




1 



(k) 



(27T) 



2</ 



/ 



z (k') i (k") i (k w ) dk'dk". 



k',k"eM M ;k'+k"+k"'=k 

Since the factor jjn is uniformly bounded and smooth for |k| 7^ conditions ()2.26|) and (|2.28j) 
are satisfied. The eigenvalues and corresponding eigenvectors of L are given explicitly: 



uu + (k) = |k| , ^ (k) = - |k| , g + (k) = 2- 1 ' 2 



I- (k) = 2- 1 / 2 



(2.101) 



Since the matrix L (k) is Hermitian, Condition 12.51 is satisfied. The singular set a consists 
of the single point k = 0. Note that conclusions of Theorem 12. 191 are applicable to equation 
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(|2.1(J(J|) and consequently to (|2.99|) . For instance, we take the initial data for (j2.1(J(J|) in the 
form ([£711 



(r, 0) = ho, V (r, Q) = h, hj = *ji (J3t - r,) e ik *< r + cc, j = 0, 1, (2.102) 



i=i 



where \l/oz ( r )j ( r ) are arbitrary Schwartz functions, cc means complex conjugate to the 
preceding terms. The points 17 are arbitrary. Note that terms corresponding to k*z can be 
written using the basis (j2.101|) as 



*0I 



e lk ^ r = [$+,/g+ + $_,,g_] e lk - r . (2.103) 



In this case all requirements of Definition 12.91 are fulfilled. The number of initial wavepackets 
for the first-order system ()2.100|) corresponding to initial data ()2.102|) ) equals Nh = 2n^ and 
there are 2Nh wavepacket centers $k^, ■& = ±. To satisfy the requirements of Condition 12. 141 
we have to take the wave vectors k^ ^0 so that 

f^^if artf^', 



\Ki' 



which provides (j2.41|) . Since 

\Ki\ A - C |Ak^| = |k^| (A - C |A|) , 

equation (|2.42j) has solutions A ^ ( and every point k^ does not satisfy Condition 12.131 This 
is the property of the very special, purely homogeneous uo (k) = |k|. Checking the second 
alternative, namely Condition 12.231 we observe that 

r-j i i i vKi v k^ 

V k WKi\ 



Hence, if 



\ukj \v\ \Ki\ 



i9k i ifk j, 

— ^ ^ —4 for I I' ot ■& ± and if k^ ^ (2.104) 



^*i\ 



\Ki' 



then Condition ^. 231 is satisfied and Superposition Theorem 12 .191 is applicable. As a corollary 
of Theorem 12.191 applied to ()2.99|) we obtain that if the initial data for ()2.99|) equal the sum 
of wavepackets, then the solution equals the sum of separate solutions plus a small remainder, 
more precisely we have the following theorem. 

Theorem 2.26 (superposition principle for wave equation) Assume that the initial data 
for i2.10(J\) to be a multi-wavepacket of the form \2.102}) and \2.4b\j holds. Then the solution 
z (r, r) to \2.1Q0i) . \2.102}) satisfy the superposition principle, namely 

z ( r > r ) = z #> i ( r ' r ) + Di ( r ' r ) ' ^ ( r ' r ) = p*.' ( r ' r ) + D2 ( r ' r ) 

■&=± 1=1 0=± 1=1 
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where z#i (r, r) 



(r, r) zs a solution of A2.1UU\) with the initial condition 
-- $^ (/3r - n) g^e ik * rr + cc, 



z#,i (r, 0) 
Pd,l (r,0) 



(2.105) 



rnt/i (r) feeing arbitrary Schwartz functions. If \2. 1 U4 ) holds, the coupling terms Di and 
D 2 satisfy the bound 



sup [HD^r) 



!D 2 (r)|U<^-^, 



0<T<T» (J 

where r* and C' s do not depend on f3,g and 5 can be taken arbitrary small. 



(2.106) 



In the following sections we introduce concepts and develop analytic tools allowing to 
prove the approximate linear superposition principle as stated in Theorems 12.151 12.191 and 
I2~25l 



3 Reduced evolution equation 

Since the properties of the evolution equations ()2.3|) and (|2.61|) are very similar, we consider 
here in detail the lattice evolution equation ()2.3|) with understanding that all the statements 
apply to the PDE (|2.61|) if we replace U with U, [— n,7i] d with R d , the function space 

L\ = L\ ^[— 7r, 7r] d j with L\ = L\ (lR d ) and so on. 

First, using the variation of constants formula we recast the modal evolution equation 
()2.3|) into the following equivalent integral form 

U (k, r) = J e Z± ? LlL(k) F (u) (k, r) dr' + e^ L(k) h (k) , r > 0. (3.1) 

Then we introduce for U (k, r) its two-time-scale representation (with respectively slow and 
fast times r and t = z ) 

U (k, r) = e-^ L{k) G (k, r) , U n , c (k, r) = u n , c (k, r) e~^" (k) , (3.2) 

where u n ^ (k, r) are the modal coefficients of u (k, r) (see (I2.18J) ): note that u n ^ (k, r) may 
depend on g, therefore (|3.2j) is just a change of variables. Consequently we obtain the 
following reduced evolution equation for u = u (k, r), r > 0, 

rap 

u (k, t)=T (u) (k, t) + h (k) , T (u) = J] ^ (m) (* m (k, r)) , (3.3) 

m=2 

jr{ m) (~mj (k? r) = £ e ^L(k)^(m) ^ ( e ^ L() a) J (k, r') dr', (3.4) 
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where the quantities are defined by (|2.21j) and (|2.22j) in terms of the susceptibilities 

The norm of the oscillatory integral jF^ m ) in ()3.4|) is estimated in terms of the norm of 

the tensor (k, k^j defined in (|2.26|) . (|2.27|) . The operator is shown to be a bounded 

one from (E) m into E, see Lemma f5. II for details The proof of this property is based on the 
following Young inequality for the convolution 

||u* v|| £l < ||u|| Li ||v|| Li . (3.5) 

For a detailed analysis of solutions of (|3.3|) we recast the equation (|3.3|) for u (k, r) using 
projections (j2.19|) as the following expanded reduced evolution equation 

oo 

< c (k, r) = E T tU (am) (k ' T) + h "<< (k) ' T - °' (3 - 6) 

m = 2 n,C 

for the modal coefficient u n ^ (k, r). In the above formula and elsewhere we use notations 



n= (n',...,n< m >), C= (C,-..,C (m) ), k= (k',...,k( m )). (3.7) 
The operators T ^ ^ are m-linear oscillatory integral operators defined by the formulas 

(Si • • ■ fi») (k, r) = f / exp (i^ (k, fc) ^ } (3.8) 

(k, £) [u! (k', n) , . . . , u m (k< m > (k, £) , 



A n,C,n,C 



where we use notations (|2.23)l . ()2.24j) . (|2.25jl . In ()3.8j) the interaction phase function (ft is 
defined by 

4> n ^ nZ (k, fc) = C^n (k) - CW (k') - ... - C (m W™> (k (m) ) , k<™> = k<™> (k, k) (3.9) 

and the susceptibilities X^^fk, are m-linear symmetric tensors (i.e. mappings from 
(C 2J ) m into C 2J ) defined for almost all k, k by the following formula 

X^ A? (k, fc) [ii! (k') , . . . , u m (k< m >)] = (3.10) 

n„, f (k) x (m) (k, fc) [n n , c (k') uj (k') , . . . , n n(m) f(m) (k< m > (k, £)) u m (k^ (k, £))" . 

For the lattice equation x*"^ - £ (k, A;^ is 27r-periodic with respect to every variable k, k', . . . , k^ m \ 
Note that operators J^i 171 ' (u m ) in ()3.3|) can be rewritten using (|3.8|) as 

We also call operators T } ^ - decorated operators. 
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Remark 3.1 The expanded reduced evolution equation \3. b]) is instrumental to the nonlinear 
analysis. Its very form, a convergent series of multilinear forms which are oscillatory inte- 
grals 113. fy) . is already a significant step in the analysis of the solution accomplishing several 
tasks: (i) it suggests a constructive representation for the solution; (ii) every term J 7 ^?.,- 
can be naturally interpreted as nonlinear interaction of the underlying linear modes; (Hi) 
the representation of as the oscillatory integral A3.8\) involving the interaction phase 

c n c an d susce Ptibilities X^ H f k, kj directly relates ^ t° the terms of the orig- 

inal evolution equation as well as to physically significant quantities. We can also add that 
since we consider g — > the interaction phase function n( -^^ (k, ttj plays the decisive role 
in the analysis of nonlinear interactions of different modes. 

The analysis of fundamental properties of the reduced evolution equation ()3.6|) . including, 
in particular, the linear modal superposition principle, involves and combines the following 
three components: (i) the linear spectral theory component in the form of the modal decom- 
position of the solution and introduction of wavepackets as elementary waves; (ii) function- 
analytic component which deals with the structure of series similar to the one in (13 .6J) and its 
dependence on the nonlinearity of the original evolution equation; (iii) asymptotic analysis of 
oscillatory integrals (|3.8|) which allows to estimate the magnitude of nonlinear interactions 
between different modes and, in particular, to show that generically different modes almost 
do not interact leading to the superposition principle. 

Sometimes it is convenient to rewrite ()3.8|) in a slightly different form. The convolution 
integral ()3.8|) according to ()2.25|) involves the following phase matching condition 

k' + ... + k (m) = k. (3.12) 

Using the following notation for the integral over the plane ()3.12|) 



/ 



/ dk' ... dk (m - 1} = (3.13) 

k',...,k( m - 1 >e[-7r,7r] (m - 1)ti ;k' + ...+k( m )=k 

/ / (k, k) 6 (k - k' - . . . - k (m) ) dk' . . . dk (m) 
in terms of a delta-function we can rewrite (|3.8|) in the form 

^ ? (u 1 ...u m )(k,r)= - * r [ e W L LRl (k,k)^\ (3.14) 



5 (k - k' - ... - k<™>) X ™tf (k, k) u 1>c (kO • • • u m>((m) (kM) dk' . . . dk^dn. 



4 Function-analytic operator series 

In this section necessary algebraic concepts required for the analysis are introduced. We study 
the reduced evolution equation ()3.3j) as a particular case of the following abstract nonlinear 
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equation in a Banach space 



u = .F(u)+x, F(u) = ^j^(x s ) (4.1) 

where the nonlinearity T (u) is an analytic operator represented by a convergent operator 
series. It is well known (see j2H]) that the solution u = Q (x) of such equation can be 
represented as a convergent series in terms of m-linear operators Q m which are constructed 
based on T: 

oo 

g(x) = g (r, x) = £ (m) ( x ™) » £ (m) ( x ™) = £ (m) x ") > where 



m=l 



X™ = X ... X = X . . . X. 



m times m times 

Using the multilinearity of we readily obtain the formula 



Q ( Xl + . . . + XJV ) = ]T £ (m) (( Xl + . . . + x^) m ) (4.2) 

m=l 

oo oo 
£ ((X^) + . . . + $ (M" 1 ) + QCI (Xi, . . . , XJV ) , 



m=l m=l 



where x = xi + . . . + x./v represents a multi-wavepacket and £?ci ( x i, • • • , x v) collects all "cross 
terms" and describes the "cross interaction" (CI) of involved wavepackets Xi,...,Xjv. We 
will find in sufficient detail the dependence of the solution operators Q m on the nonlinearity T 
and prepare a basis for the consequent estimation of nonlinear interactions between different 
modes and wavepackets. Then combining the facts about the structure of the solution 
operators with asymptotic estimates of relevant oscillatory integrals we show that for a 
multi-wavepacket x = xi + . . . + x^ the cross-interaction term satisfies the following estimate 

H0CI (xi, . . . , XiV ) || = O (/?) + O (q |ln /3| , 0,q - 0, 

implying the modal superposition principle. 



4.1 Multilinear forms and polynomial operators 

The analysis of nonlinear equations of the form (|3.3|) requires the use of appropriate Banach 
spaces of time dependent fields, as well as multilinear and analytic functions in those spaces. 
It also uses an appropriate version of the implicit function theorem. For the reader's conve- 
nience we collect in this section the known concepts and statements on the above-mentioned 
subjects needed for our analysis. We in this section consider functional-analytic operators 
which are defined in a ball in a Banach space X with the norm ||x|| x . In our treatment of 
the analytic functions in infinitely-dimensional Banach spaces we follow to [23 Section 3], 

EH- 
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Definition 4.1 (polylinear operator) Suppose that x 1; x 2 , . . . , x n are vectors in a Banach 
space X. Let a function T^ n ' (x), x = (x 1; . . . ,x„), take values in X and be defined for all 
x G X n . Such a function is called a n-linear operator if it is linear in each variable, and 
it is said to be bounded if its following norm is finite 

||^ {n) ||= sup ||^(x 1 x 2 ...x n )|| x < oo. (4.3) 

ll x l \\x = --- = ll x ™ II X = l 

Definition 4.2 (polynomial) A function P (x) from X to X defined for all x G X is called 
a polynomial in x of degree n if for all a, h G X and all complex a 

n 

P(a + ah) = ^P y (a,h) a\ 

where P v (a, h) G X are independent of a . The degree of P n is exactly n if P n (a, h) is not 
identically zero. A polynomial T (x) is a homogeneous polynomial of a degree n if for all 
CG C 

T (cx) = c n T (x) . 

Then n is called also the homogeneity index o/jF(x). A homogeneous polynomial J 7 is called 
bounded if its norm 

OT.= ^p {||^(x)|| x } (4.4) 

IMIx =1 

is finite. For a given n-linear operator (x) = T^ n * (xix 2 . . . x n ) we denote by (x n ) 
a homogeneous of degree n polynomial from X to X: 

jF(")( x ") = ^ (n) (x...x). (4.5) 

Note the norm definitions ()4.3|) . ()4.4|) and ()4.5|) readily imply 

||jrW||^ < ||jr(«)|| , ( 4 6 ) 

Definition 4.3 (analyticity class 1) Let a function T be defined by the following conver- 
gent series 

oo 

F{x) = J2^ m) (x m ) for ||x|| x < R^, (4.7) 

m=2 

where (x m ) ; m = 2,3, ... is a sequence of bounded m-homogenious polynomials satisfying 

\\F {m) \l<C^R;™, m = 2,3,..., (4.8) 
Then we say that T (x) belongs to the analyticity class A* (C*jf, R^) and write T G A* (C^, R*?) . 

Notice that for ||x|| x < R^ we have 

00 ll-vll n ° R~ n ° 

n=2 1 \\ X Wx n *F 

implying, in particular, the convergence of the series ()4.7j) . 
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Definition 4.4 (analyticity class 2) If T^ m ' (x), m = 2, 3, . . is a sequence of bounded 
m-linear operators from X m to X and 

||^ (m) || <C T R~ m , m = 2,3...., (4.10) 

we say that a function T defined by the series \4- f or \\ x \\x < belongs to the analyticity 
class A (Cy, Rjt) and write T G A (Cjr, R^) . 

In this paper we will use operators from the classes A (Cjr, Rj?) based on multilinear 
operators. 

Note that evidently A (CV, Rp) C A* (Cy, Pp). One can construct a polynomial based 
on a multilinear operator according to the formula (|4.5jl . Conversely, the construction of a 
multilinear operator, called polar form, based on a given homogeneous polynomial is described 
by the following statement, [2H Section 1.1, 1.3], [221 Section 26.2]. 

Proposition 4.5 (polar form) For any homogeneous polynomial (x) of degree n there 
is a unique symmetric n-linear operator (X1X2 . . . x„,) ; called the polar form of P n (x), 
such that P^ (x) = P^ (x . . . x). It is defined by the following polarization formula: 

PW (x lX2 . . .x„) = JL Y: P W (Ee^) . (4.11) 

' e,-=±i Vj=i / 

In addition to that, the following estimate holds: 



< 



Pn 



<^\\Pn\l<e n \\P {n) L- (4-12) 



Since by Definition 14.41 functions from A (C, R) have zero of the second order at zero, 
their Lipschitz constant is small in a vicinity of zero. More exactly, the following statement 
holds. 

Lemma 4.6 (Lipschitz estimate) // T E A(Cp, R?) then 

||^(x)-^(y)|| <CyCJx-y||(||x|| + ||y||) for ||x|| , ||y|| < Rl T < R T , (4.13) 
where C > depends on R! T and Rjr. 



4.2 Implicit Function Theorem and expansion of operators into 
composition monomials 

Here we provide a version of the Implicit Function Theorem, first we formulate classical 
implicit function theorem for equations u = T (u) + x with analytic function JF and then we 
present a refined implicit function theorem. The refined implicit function theorem we prove 
here produces expansion of the solution u into a sum of terms which are multi-linear not only 



36 



with respect to x but also with respect to T . The formulation of the theorem and the proof 
involve convenient labeling of the terms of the expansion (called composition monomials), 
and we use properly introduced trees to this end. The explicit expansion produced by the 
refined implicit function theorem is required to be able to take into account rather subtle 
mechanisms which lead to the superposition principle. 

Let us consider the abstract nonlinear equation (J4.1)) and its solution u = u (x) for small 
||x|| when the nonlinear operator JF belongs to the class A(Cf, RjA. We seek the solution u 
in the following form 

oo 

u = g {T, x) = £ M ( x ™) for sufficiently small ||x|| , (4.14) 

m=l 

and we call Q the solution operator for (J4.1j) . It readily follows from (J4.1j) that 

0(.F,x) = x + .F(S(jF,x)) (4.15) 

and 

oo oo / / oo \ s \ 

j2 Q [m) (x m ) = x + f {3) Yl g{m) (* w ) • ( 4 - 16 ) 

m=l s=2 \ \m=l / / 

^From the above equation we can deduce recurrent formulas for multilinear operators g( m \ 
Indeed for m = 1 the linear term is the identity operator 

g(V ( x ) = JFW ( x ) = x. (4.17) 
For m > 2 we write the following recurrent formula 

m 

( Xl ...x m ) = £ Yl ^ ( g{h) ( x i • • • x u) • • • Q [ta) (x m _, s+ i • • • x m )) . (4.18) 

s=2 i\+...+i a =m 

By the construction, if multilinear operators are defined by (|4.18|) . then (|4.16|) is satisfied. 
Namely, expanding right-hand side of (J4.16)) using multi-linearity of JF^ we obtain a sum 
of expressions as in right-hand side of (|4.18|) . and since ()4.18|) holds, terms in the left-hand 
side of ()4.16j) with given homogeneity index p cancel with the terms in the right-hand side 
with the same homogeneity. Note that in ()4. 18|) we do not assume that the operators 
and Qv) are symmetrized and the order of variables is important; we prefer to treat 
and g( m ^ as multilinear operators of s and m variables respectively. Though, when we apply 
constructed to solve (|4.1|) . we set xi = . . . = x m . 

The following implicit function theorem holds (see [3] and Theorem 14.251 below with a 
similar proof). 

Theorem 4.7 (Implicit Function Theorem) Let T G A[Cjr^R :F ). Then there exists a 
solution u = x + g (JF, x) of the equation \4-l]) u = x + T (u) ; given by the solution operator 
g E A (Cg, Rg), where we can take 

R 2 R 2 

Cg = 2{c/ +Rr y Rg = mc^rVV (419) 
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the series \4-14\) converges for ||x|| x < Rg. The multilinear operators Q^ m > (x) satisfy the 
recursive relations j^. \4-l$ - 

Note that uniqueness of the solution and continuous dependence on parameters follows 
from Lemma I4~fi1 and from a standard observation which we formulate in the following remark. 

Remark 4.8 Ifii\, U2 are two solutions of the equation \4-l\) with x = xi, X2 respectively and 
|| u i|| ? || u 2|| < R, and ^"(u) is Lipschitz continuous for ||u|| < R with a Lipschitz constant 
q < 1 then \\u% — u 2 || < (1 — q)^ 1 ||hi — h 2 || . If U!,u 2 are two solutions of the equation 
( |^.i| ) with T = T§ and T = J-q + T\ respectively, ||ui|| , ||u 2 || < R, and T {vl) is Lipschitz 
continuous for ||u|| < R with a Lipschitz constant q < 1 and T\ (u) < e when ||u|| < R then 
||ui - u 2 || < e (1 - q)' 1 . 

Observe that every term in (J4.18j) . in turn, can be recast as a sum (|4.18|) with m 
replaced by ii < m. Evidently applying the recurrent representation (J4.18)) and multilinearity 
of ;F( S ) we can get a formula for Q^ m ' as a sum of terms involving exclusively (i) the symbols 
JFM, (ii) variables Xj and (hi) parentheses. We will refer to the terms of such a formula 
as composition monomials. To be precise we give below a formal recursive definition of 
composition monomials. The monomials are expressions which involve variables uj, j = 
1,2,..., and m-linear operators JF*" 1 ), m = 2,3..., and are constructed by induction as 
follows. 

Definition 4.9 (composition monomials) Let {-T 7 ^}^ be a sequence of s -linear oper- 
ators which act on variables Uj, j = 1,2,.... A composition monomial M of rank is 
the identity operator, namely M (uj) = Uj, and its homogeneity index is 1. A composition 
monomial M of a non-zero rank r > 1 has the form 

M (u l0 . . . iiij = ^ (s) (Mi (u J0 . . . UiJ . . .M s (u ls _ 1+1 . . . ii,.)) , (4.20) 

where M x (u io . ..u h ), M 2 (u il+1 . . . u i2 ),..., M s (u is __ 1+1 . . .u is ), with 1 < i < % x < . . . < i s , 
are composition monomials of ranks not exceeding r — 1 ( submonomials ) and at least one of 
the rank r — 1, the homogeneity index of Mj equals ij — ij-i- For a composition monomial 
M the operator JF^ in its representation \4-%$ is called its root operator. The index of 
homogeneity of M defined by fl^.£fl| ) equals i m — io + 1. We call the labeling of the arguments 
of a composition monomial M defined by M-^Uj ) by consecutive integers standard labeling if 
io = 1- 

If the monomials M l5 .,M s have the respective homogeneity indexes v {MA then we 
readily get that the homogeneity index of the monomial M satisfies the identity 

v (M) = v (Mi) + ... + V (M s ) . (4.21) 

Using the formula (J4.20)) inductively we find that any composition monomial M is given by a 
formula which involves symbols from the set {.F (s ' l }^ 2 , arguments Uj and parentheses, and 
if s-linear operators are substituted as J 7 ^ we obtain the terms contained in the expansion 

of g {m \ 
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Definition 4.10 (incidence number) The total number of symbols involved in M is 
called the incidence number for M. 

For instance, the expression of the form 

M = JF (4) (u 1 u 2 u 3 ^ {3) (u 4 ^ {2) (u 5 u 6 ) (u 7 u 8 u 9 )) ) (4.22) 

is an example of a composition monomial M of rank 3, incidence number 4 and homogeneity 
index 9. It has three submonomials. Namely, the first one is J 7 ^ (u 4 JF( 2 ) (u 5 u 6 ) JF® (u 7 u 8 u 9 )) 
of rank 2 and incidence number 3. The second submonomial (u 5 u 6 ) has rank 1 and in- 
cidence number 1, and the third one is 

^(3) 

(u 7 u 8 u 9 ) of rank 1 and incidence number 1. 
When analyzing the structure of composition monomials we use basic concepts and no- 
tation from the graph theory, namely, nodes, trees and subtrees. 

Definition 4.11 (nodes, tree, subtree) A (finite) directed graph T consists of nodes 
Ni e Nt where is the set (finite) of nodes ofT and a set of edges NiNj e x N?- An 
edge N^Nj connects Ni with Nj, it is an outcoming edge of Ni and an incoming edge of Nj. A 
tree (more precisely a rooted tree, we only consider rooted trees) is a directed connected graph 
which is cycle-free and has a selected root node, that is a node N* which has no incoming 
edges. If a node N has an outcoming edge NNj the node Nj is called a child node of N; if a 
node N has an incoming edge NjN the node Nj is called the parent node of N. We denote 
the parent node of N by p (N) . If a node does not have children it is called an end node ( or 
a leaf). For every node N we denote by /i (N) the number of child nodes of the node N. If 
a path connects two nodes we call the number of edges in the path its length. We denote 
by I (N) the length of a path which connects iV* with N. Every node N of the tree T can be 
taken as a root node of a subtree which involves all descendent nodes of N and connecting 
edges; we denote this maximal subtree T' (N) . Since we consider only maximal subtrees we 
simply call them subtrees. We call by the rank of a tree the maximal length of a path from 
its root node to an end node and denote it by r (T) . We call by the rank of a node iV of the 
tree T the rank of the subtree T' (N) . 

Definition 4.12 (tree incidence number and homogeneity index) For a treeT we call 
the number of non-end nodes incidence number i — i (T). We denote the number of end nodes 
of the tree by v (T) and call it homogeneity index. 

Elementary properties of trees. Since a tree does not have cycles, the path connecting 
two nodes on a tree is unique. The root node iV* does not have a parent node, and since it 
is connected with every other node, every non-root node has a parent node. The end nodes 
have zero rank. The only node with rank r (T) is the root node. The total number of nodes 
of a tree T equals m{T)+i (T). 

Definition 4.13 (ordered tree) A tree is called an ordered tree if for every node N all 
child nodes of N are labeled by consecutive positive integers (which may start not from 1). 
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Hence, for any node N' ^ JV* there is the parent node N = p (TV 7 ) and the order number (label) 
o(N'), %x < o(N') < %x + (J>(N) — 1. Two trees are equal if there is one-to-one mapping G 
between the nodes which preserves edges, maps the root node into the root node and preserves 
the order of children of every node up to a shift: if (JV) = N and p (JVi) = p (JV2) = N 
theno{N l )-o{N 2 ) = o (9 (JV^) - o (6 (JV 2 )). 

Since we use in this paper only ordered trees we simply call them trees. 

Standard node labeling and ordering. We use the following way of labeling and order- 
ing of end nodes of a given ordered tree T. Let f be the rank of T. For any end node N we 
take the unique path N*Ni . . . N^ N )^iN of length I (JV) < f connecting it to the root. Since 
the tree is ordered, every node Nj in the path has an order number o (Nj). These order num- 
bers form a word w (N) of length I (JV). If / (JV) < f we complete w (JV) to the length f adding 
several symbols 00 and assuming that 00 > n for n — 1,2, . . .. After that we order words 
w (JV) in the lexicographic order. We obtain the ordered list w 1 (Ni) , . . . ,w u ( T ) (N u (t)). We 
take this ordering and labeling of the end nodes N ± , . . . ,N v (j>) as a standard ordering and 
denote by Oq (N) the consecutive number with respect to this labeling: j = oq {Nj). To label 
the nodes with rank r we delete all the nodes of rank less than r together with the incoming 
edges and nodes of rank r become end nodes. We apply to them the described labeling 
and denote the indexes obtained by o r (N). Hence, every node N of the tree T has two 
integer numbers assigned: r (N) and o r ( N ) (N). We introduce the standard labeling of all 
nodes of T by applying the lexicographic ordering to pairs (r (JV) , o r (N) (iV)) , and denote the 
corresponding number o (N), 1 < o (N) < m (T) + % (T). 

The following statement follows straightforwardly from the definition of the standard 
ordering. 

Proposition 4.14 If a tree T has a subtree T' and the standard labeling of end nodes is 
used, then all the end nodes of the subtree T' fill an interval j\ < o (JV) < j 2 for some j\ 
and J2 ■ 

Theorem 4.15 Let T2 be the set of ordered trees such that each node of a tree which is 
not an end node has at least two children nodes. The set of composition monomials based 
on {jF( s ), s = 2, 3, . . .} is in one-to-one correspondence with the set T 2 . The correspondence 
has the following properties. The monomials of rank r correspond to trees of rank r. The 
root node of the tree T corresponds to the root operator of the composition monomial. The 
end nodes correspond to variables Uj, j = 1, . . . , v (T). The standard labeling of end nodes 
coincides with the consecutive labeling of the variables \\j of monomial from left to right. The 
homogeneity index of a monomial equals the homogeneity index of the corresponding tree. 
The incidence number of a monomial equals the incidence number of a tree, and the rank of 
a monomial equals the rank of a tree. 

Proof. For a given {jF^* 1 } the set of monomials with rank r is finite, the set of trees 
with rank r is finite too. Therefore, to prove one-to-one correspondence of the two sets it is 
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sufficient to construct two one-to-one mappings from the first set into the second and from 
the second into the first. First of all, using the induction with respect to r we construct for 
every monomial the corresponding tree. Let r = 0. A monomial of rank has the form Ui, 
and it corresponds to a tree involving one node. The tree has no edges and the node is the 
both the root and the end node; its incidence number is zero and homogeneity power is one. 
Assume now that we have defined a tree for any monomial of rank not greater than r — 1. A 
monomial of rank r has the form JF^" 1 ) (Mi . . . M m ) where monomials Mi . . . M m have rank 
not greater than r — 1 . Every monomial Mi . . . M m corresponds to an ordered tree T% , . . . , T m 
with the root nodes iV*i, . . . N* m . We form the tree T as a union of the nodes of Ti, . . . , T m 
and add one more node iV* which corresponds to the root operator JF*™) and it becomes 
the root node of T. We take the union of edges from T 1; . . . ,T m and add m more edges 
connecting N* with the nodes iV*i, . . . N* m , the order of the nodes corresponds to ordering of 
Mi . . . M m from left to right. The first mapping is constructed. 

Now let us define for every ordered tree T the corresponding monomial M (JF, T). If we 
have a tree T of rank zero we set M (JF, T) = u, and j — 1 if we use the standard labeling. 
Now we do induction step from r — 1 to r. If we have a tree of rank r we take the root node 
iV* and its children N*i, . . . , N* s , s — fi (iV*). The subtrees T' (N*i ),..., T' (N* s ) have rank not 
greater than r — 1 and the monomials M (JF, T' (-/V*i)) M (J 7 , T' (N* s )) are defined accord- 
ing to induction assumption, let m (T f (JV*i )),..., m (T' (N* s )) be their homogeneity indices. 
We set m (T) = m (T" (iV*i)) +...+m (T 1 (iV* s )). We denote the variables of every monomial 
M (J 7 , T' (iV*j)) by Uj i, . . . , \ij im (T'(N t ■)) counting from left to right, and then labeling all the 
variables ixy using the lexicographic ordering of pairs j, I we obtain variables u 1; . . . , u m (r) 
and monomials 

M (J 7 , V (N^)) (ui, . . . , u m (T/(jv» 1 )) ) , M (J 7 , T' (iV* 2 )) ( 

etc., where m 3 - = m (T' (N*j)). After that we set 

M(F,T) (ui, . . . ,u m(r) ) = 
(M (^, T' (N*i)) (ui, . . . , ^/(jv.o)) , • • • , M (JF, T' (JV„)) (u m( T)-m._ 1 +i, • • • , u m( r))) . 

Note that the homogeneity index for the monomial M equals the sum of the indices for 
submonomials Mi . . . M m , the homogeneity index for the tree T equals the sum of the indices 
for subtrees Ti, . . . , T m , this implies their equality by induction. The incidence number for 
the monomial M equals the sum of the numbers for submonomials Mi . . . M m plus one; the 
incidence number for the tree T equals the sum of the numbers for submonomials Ti , . . . , T m 
plus one. Therefore, these quantities for monomials and trees are equal by induction. 
Induction is completed. Therefore we constructed the two mappings, one can easily check 
that they are one-to-one and have all required properties. ■ 

Definition 4.16 (monomial to a tree) For a tree T E T 2 we denote by M(jF, T) the 
monomial which is constructed in Theorem \4- 73| 
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Conclusion 4.17 The above construction shows that the structure of every composition 
monomial is completely described by an (ordered) tree T with nodes iVj corresponding to 
the operators J r ^ nii \ At such a node Ni (i) the number rrii of outcoming edges equals the 
homogeneity index of J 7 ^' ; (ii) the outcoming edges are in one-to-one correspondence with 
the arguments of J r ( rrii \ and the ordering of the child nodes coincides with the ordering of 
arguments of J 7 ^ 71 *) from left to right. The value of mi may be different for different nodes. 
A node corresponding to J 7 ^ is connected by edges with m child nodes corresponding to the 
arguments of J 7 ^ . Every node N of the tree T can be taken as a root node of a subtree 
T' (N) which correspond to a submonomial M (J 7 ,T' (N)) . Conversely, every submonomial 
of M {? ', T) equals M (J-,T' (N)) for some mode N. If m > 1 the submonomial has a 
nonzero rank. The number of non-end nodes equals to the number of symbols J 7 ^ used in 
T -represenation of the monomial which is the incidence number of the monomial. The total 
number of end nodes of an m-homogenious operator equals to m = v (T). The rank of a node 
N equals the rank of the corresponding submonomial M (J 7 , T' (N)). The arguments ui, . . . u s 
of a monomial correspond to the end nodes of the tree. The standard labeling of nodes of 
T agrees with the standard labeling (from left to right) of the arguments of the composition 
monomial M (J 7 , T). The number of end nodes of the tree T equals the homogeneity index of 
corresponding monomial. If the root mode of the tree T of a monomial M has [i (iV*) = m 
edges which are connected to child nodes N%, . . . N m then there is a node JF( m J') ; j = 1, . . . , n 
at the end of every edge such that M has the form 

T {m) ^(JVO) (...),..., jrMN m )) _ ^ _ ( 423 ) 

Example 4.18 The tree corresponding to J 7 ^ (u 1 u 2 JF (111U2U3)) has two nodes of non-zero 
rank, the root node of rank 2, one non-end node of rank 1 and five end nodes of rank 0. 
Another example, the monomial \4-2^j) has the root node corresponding to J 7 ^' , four edges 
lead respectively to nodes corresponding to the end nodes with ui, U2 ; U3 and to the non-end 
node with J 7 ^ , see Fig. 2. 

Remark 4.19 Since all operators in the set {-T 7 ^}^ in \4-l£fy have the homogeneity index 
at least two, the trees of monomials generated by recurrent relations \4.lty) have a special 
property: every non-end mode has at least two children. 

Sometimes it is convenient to use monomials involving several types of operators. To 
describe such a situation we introduce for a given tree a decorated monomial. 

Definition 4.20 (decorated monomial of a tree) Assume that we have several formal 
series . . . , J 7 ]] where J 7 ^ is represented by a formal series T\ = J 7 ^, % = 1, . . . , I. 
We call the set {J 7 } = {J 7 ^, j — 1, . . . , S} the operator alphabet, and j is called the deco- 
ration index. We consider a function T (N), N G T, defined on the nodes of the tree T and 
taking values in the set {1, . . . , /} of the decoration indices, and call such a function a deco- 
ration function on the tree T. Then for a decoration function F (N) we define the decorated 
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F (4) (uiu 2 u 3 F (3, (U4 F (2) (u 5 u 6 )F (3, (u7U S u 9 ))) 



Figure 2: In this picture a tree corresponding to a monomial is drawn. 

monomial M ({JF} , T, T) of the tree T by picking operators J-^ with j defined by T. For 

every node N the homogeneity index m — // (N) of the operator equals to the number 
of children of N and j is defined by T, namely J-j, j = T (N). 

Hence, a decorated monomial M {{J 7 } ,T,T) has instead of ()4.23|) the following form 

-r-M ( -r-MNi)) / \ T (n(N m )) i \\ ( A 0A \ 

• r r(Af) ^* r r(Af i ) \- ■ ■> ' • • • ' • r r(Ar m ) v • •)) ■ \^-^> 

When are multilinear operators, a monomial M({JF},T, T) is also a multilinear op- 

erator, its homogeneity index m equals v (T) and we denote its arguments by (xi . . . x m ) . 
Respectively, if xi . . . x„ are arguments of a monomial M {{J 7 } , T, T) and we use the stan- 
dard labeling of the nodes then according to Proposition I4.14| a submonomial M {{J 7 } , T, V) 
has arguments x^/), . . . , x*c(t')+i>(t')-i which are labeled constructively 

Now we would like to describe elementary properties of composition monomials and the 
related trees. Note that for every iV e T a composition monomial is a linear function of 
operator ^(ao- Consequently, the concept of the decorated composition monomial can be 
naturally extended to monomials associated with the following family of operators 

{F} = {F:r=c 1 r 1 + ... + cfa c t eC}. 

For a given tree T the submonomial M ({J 7 } , T, T) is represented as a function on the tree 
T with values in {J 7 }, this is an i-linear function of T where % is the incidence number of T. 
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There are elementary relations between the incidence number i (T), the rank r (T), the 
number of edges of a tree T which do not end at an end node e° (T) and the homogeneity index 
m of a tree T, and corresponding monomial M ({J 7 } , T, T). For example, e° (T) = i (T) — 1. 
Some useful relations expressed by inequalities are given in the following lemma. 

Lemma 4.21 Let us consider trees T for which every non-end node has at least two children, 
^ (N) > 2 for all N e T. Let for any i the number m (i) be the minimum number of the end 
nodes v (T) for all trees T with given incidence number i. Then 

m(i)>i + l. (4.25) 

Similarly for any given r let m (r) be the minimum number of end nodes with given rank r. 
Then 

m(r)>r + l. (4.26) 

Let e° (T) be the number of edges of a tree T which do not end at end nodes. For any given 
e let m (e) be the minimum number of end nodes with e° (T) = e. Then 

m (e°) > e° + 1. (4.27) 

Proof. For i — 1 (J4.25|) is true. Let the statement be true for i = i . Let T be a tree with 
the minimum number of end nodes m (io) = m. We delete one of the end nodes together 
with the edge leading to it from its parent obtaining a tree with m (io) — 1 end node. If the 
tree remains in the same class, then m (io) is reduced by one contradicting the minimality. 
Hence, the deletion of the edge created a node with only one child. Such a node can be 
replaced by an edge leading from its parent to its child and reducing the incidence number 
by one. Using the induction assumption we get 

m(i ) - 1 > m(z - 1) > (i - 1) + 1 (4.28) 

that completes the induction and proves (|4.25|) for all i. Similar induction proves (|4.26|) . For 
r = 1 ()4.26|) is true. Let T be a tree with the minimum number of end nodes m (ro) = m. 
As above, by deleting an end node and using the minimality we reduce the tree T to a tree 
T' with a smaller rank. Since only one non-end node is eliminated, the rank of T' is ro — 1 
and we get (|4.26|) . Inequality (J4.27|) holds for e = since m (0) > 2. Let T be a tree with the 
minimum number of end nodes m (eo) = m. We again delete one of the end nodes together 
with the edge joining it to its parent and obtain a tree with m (eo) — 1 end nodes and the 
same number of edges which do not end at an end node. The minimality implies that the 
parent node has only one another child and removing it we get either eo or eo — 1 edges which 
do not go to end nodes. We use the induction as in (|4.28j) obtaining (|4.27|) . ■ 

Monomial expansion in the Implicit Function Theorem If operators Q m (x x . . . x m ) 

are determined by the recurrent formulas (J4.18)) it is obvious that every Q m can be represented 
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in terms of T = {.F^} using the recurrence and multilinearity of T^ 8 ' . More precisely the 
following representation holds 

g {m) (T, x x . . . x m ) = °t M (F, T ) ( x i • • • O , (4.29) 

TeTm 

where (i) M (JF, T) G T 2 is a composition monomial corresponding to a tree T and T m C T 2 
stands for the set of trees with m end nodes; (ii) the integer-valued multiplicity coefficient 
ct > counts the multiplicity of the related monomial M(jF, T) in the expansion of (j4.18|) : 
for some trees T its multiplicity coefficient ct may be zero. The expansion ()4.29|) is obtained 
by an inductive process with respect to m since (J4.18|) expresses Q m in terms of Q^ 3 ' with 
2 < ij < m. Notice that for a given operator T = {.F^-*} the monomial M (.F, T) considered 
as an operator can be the same for different T, the monomials and the multiplicity coefficients 
are determined purely algebraically. 

Remark 4.22 The expression for as a linear combination of composition mono- 

mials M (JF, T), in particular the multiplicity coefficients ct, does not depend on a specific 
form of the operator T . It is the same for a solution z = x + Q (JF, x) of the general func- 
tional equation J^.i| ) and for an elementary algebraic equation u = T (u) + x with u, x G C 
and with a scalar analytic function T of one complex variable. 

If all T { ™ ] are bounded multilinear operators then a decorated monomial M (JF, T, T) is 
also a bounded multilinear operator as it follows from the following statement. 

Lemma 4.23 Let M ({J 7 } , T, T) be a decorated monomial of the homogeneity index v (T) = 
m and all J 7 ^ be bounded operators from E s into E for a Banach space E. Then the following 
estimate holds 

m 

iiM({^},r,r)(x 1 ...x m )n B < n n inu- ( 4 - 3 °) 

NaT ,r(N)>0 j=l 

Proof. Notice that 

\\^ {M 1 . . . M m )\\ E < H^ll HMJ^ . . . \\M m \\ E (4.31) 

where Mj are submonomials. Applying the above inequality repeatedly we obtain (j4.30J) . ■ 
The next statement provides a bound for the norm of a decorated monomial which involves 
factor the norm of a submonomial. 

Lemma 4.24 . Let M ({J 7 } ,T,T) be a decorated monomial evaluated at xi...x m . Let 
all be bounded operators from E s into Banach space E. Then for every evaluated 
submonomial M ({J 7 } , T' (Nq) , T) we have an estimate 

\\M ({J 7 } , T, T) ( Xl . . . x m ) \\ E < \\M ({J 7 } , T' (N ) , T) ( Xje , . . . , ** +v{t > { n))-i) \\ e (4.32) 

NeT\T'(N ),r{N)>Q i<x j>x+v(T>(N )) 

where x^ , . . . , x^ +u (t>(n))-i are the arguments of the submonomial M ({J 7 } , T' (Nq) , T). 
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Proof. The proof uses the induction with respect to the length I (No). For I (N ) = 
the statement is obvious. Assuming that the statement is true for / (N) < Iq we consider 
the case when I (N ) = l . Notice that 



Oi(Jv.)) 



(Mi . . . M KN) ) 



< 



\MA 



E ■ ■ ■ 



M 



where Mj = M ({J 7 } ,T' (N*j) ,T), N*j are child nodes of N*. One of the submonomials 
Mi . . . M^n) contains M ({JF} , T (N ) , T) as a submonomial, and let it be M ({^} , T (N* jo ) , T). 
The length of the path from iVo to N*j is less than l Q and we can use the induction hypothesis 
to estimate the norm of M ({J 7 } , T' (N*j ) , T). The norms of Mj with j ^ jo are estimated 
using (|4.30jl . The labels of the arguments of the submonomial fill an interval according to 
Proposition 14.141 ■ 

The following theorem gives a needed refinement of the Implicit Function Theorem 14.71 

Theorem 4.25 (refined Implicit Function Theorem) Let T E A(C?, Rr). Let Q E 
A(Cg,Rg) be the analytic solution operator constructed in Theorem \4-7\ which solves \4-l\j - 
Then the expansion of Q (J-, x) into composition monomials 



s (.f, x) = E E C ? M P> T ) ( x " 

m=l TGT m 

converges for ||x|| < Rg, and the following estimates hold 
E c T \\M(F,T)(x m )\\<C g Rg m \\x\ 



(4.33) 



m 



(4.34) 



E E c T \\M(F,T)(x m )\\<Cg- 

m=2 T&T m 

where Cg and Rg depend only on and Rjr and satisfy 



\ x \\x 



\ x \\x Rg 



C c 



R% 



JQ 2(C r + R r y Rg A(C r + R r y 
The multiplicity coefficients ct > satisfy the inequality 



R 2 T 



E c ^ ¥ m 



(4.35) 



Ter n 



The proof of this statement is given in Appendix B. 
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4.3 Decorated expansions 

In this section we develop a formalism for treating linear operators with several invariant 
subspaces which span the entire space as, for example, in the case of projections (|2.19j) . 
The decomposition into related invariant subspaces is very important for the analysis. The 
general setting is as follows. Suppose that a Banach space E has several projection operators 
IT A , A G A, where A is a finite set of indices, we call this set decoration set. We assume that 
the sum of the projections equals the identical operator, i.e. 

n A = Id, where Id is the identity operator, (4.36) 

AeA 

and 

n A n A = n A , n v n A = o if x' ^ x, a', a e a. (4.37) 

We call such projections decoration projections. For example, let us look at projections 
U n ^ (k), n = 1, . . . , J, ( = ± defined by (|2.19j) . These projections define bounded operators 
H n ,( acting on (i) functions of k in the space L\\ (ii) functions of k, r in the space E = 
C ([0, t*] , Li). In another example based on (j2.19|) we fix n and define 

n c (k) = n„ 0iC (k) , c = ±, (k) = £ n n , c (k) . (4.38) 

Using (|4.36j) we expand vectors x G E as follows 

AeA AeA 

We also use notation 

jf) = UxFW (4.40) 

Often in applications the number of elements in A is either 2 or 3. In the case when A has 
three elements we set 

a = {+, -, 00} , n + + n_ + iioo = id, (4.4i) 

and 

x = x+ + x_ + Xoo , T (x) = F + (x) + T- (x) + (x) . (4.42) 

Using the decomposition ()4.36|) we introduce for m-linear operators (xi . . . x n ) the cor- 
responding decorated operators T[ n ) as follows: 

J** ( Xl . . . x„) = U X T^ (n c / Xl . . . n c( „)X„) = T { x n) (n c , Xl . . . U, w x n ) , (4.43) 

where ( is defined in (13. 7j) . Obviously, we have 

j&> ( Xl . . . Xn ) = Yl ( x i • • • x «) ■ ( 4 - 44 ) 

AeA, ceA" 

An example of expansion (J4.44j) is given by (|3.11j) . 



47 



4.4 Decorated composition monomials 

We assume that operators J 7 ^ 71 ' act in the space allowing a decomposition into three com- 
ponents as in (|4.41|) . Let M(T,T) be a composition monomial of the homogeneity index 
m, and assume that the corresponding tree T has the incidence number i, the rank r, and e 
edges. Suppose also that every operator J 7 ^ is expanded into a sum of decorated operators 
as in (J4.44j) , (J4.43j) . Using the linearity of M (J 7 , T) with respect to operators J 7 ^ we get 



M (J 7 , T) = J 7 ^ (^ (mi) (...)... -F (mn) (• • •)) = (4-45) 

r (mi) ^ ^ jr(m n ^ ^ 

^31 'Ol Aj'n >f jn 



^ \ A-,*, A 



AeA, AeA'- 1 , Cj, j=l,...,e 

where submonomials J 7 ^ 1 -) (. . .),..., J r< - m ™- ) (...) have ranks not exceeding r — 1. We expanded 
repeatedly the expression in the left-hand side of ()4.45|) as long as submonomials of non-zero 
rank were present resulting in an expansion involving only decorated operators J-[ n ). 

Remark 4.26 Note that 



(j 7 ^ (...)••• J**? (• • •)) = ^ ( n Al ^ m i } (...)... n A „^ (. . .)) (4.46) 



Smce projections satisfy the identities \4- 31§ if a vector (f = k', . . . , <^ n ^ J and indices 
Ai, . . . , A n are given, then we have the identity 

^{n) f f W _ _ _ ^(m ) \ = Q when ^ ^ some . (4 47) 

A,C V Ai,Ci An,C„/ 

Hence, for non-zero terms in the expansion (R~73L) z/ indices Ai, . . . , A n /or jF| m ^,... ; jF| m ".' ) 

Ai,Ci A n ,C„ 

are azwen iae vector ( in J 7 ^ is determined by them 

( {l) = \ t , i = l,...,n. (4.48) 
Note that according to (J4.47|) and (J4.48)) we have 

-pH jjK) (...) ... jr(»n») (. . .)) = jr(n) fjjnx) (...)... jr<™™) (. . .) \ . (4.49) 

According to ()4.45|) , ()4.49|) for every tree T of the homogeneity index m and the incidence 
number i, we get an expansion into a sum of monomials of the form 

M {?, T, A, () (xix 2 . . . x m ) = M {{J 7 } , r, T) ( Xl x 2 . . . x m ) , (4.50) 
{^} = {^:A6A, CGA", n = 2,3,...}. 
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Namely, if a monomial M (J 7 , T) has at a node N operator J^( m ( N )) then M ({J 7 } , T, T) at 
this node has operator ^F^r^P ■ We call a composition monomial of the form (J4.5(J|) . where 
f!4.48|) is assumed, a decorated composition monomial. Using the standard labeling of nodes, 
for a given function T on the tree T with values in A we find the vectors A G A 1 , ( e A" 1 , 
with i being the incidence number of the tree T, and using (|4.48|) we rewrite (|4.45|) in the 
form 

M(^,T)( Xl x 2 ...x m ) = Yl M(T,T,\,t) (x lX2 ...x m ). (4.51) 

AeA 1 , <feA m 

where ( is determined by values of T on the end nodes. The sum (J4.51|) contains at most 
2«+m non _ zero terms, where 3 is the number of elements in A. Combining (J4.51|) with (J4.33j) 
we obtain 

g(m) (xm) = J2 J2 °tM (>, T, A, C) (x m ) . (4.52) 

TeT m \£Ai(T) t ^ eA m 

5 Expansions of solutions for oscillatory integral equa- 
tion 

In this section we apply general concepts introduced in previous sections to oscillatory in- 
tegrals involving operators T as in ()3.3|) . ()3.4|) . Based on projections n n ^ (k) in ()2.19|) for 
given n = n we define as in ()4.38|) decoration projections in L\ which satisfy ()4.41|) : 

n c u (k) = n no , f (k) u (k) , c = ±, Hoc = E n ™.c- (5-i) 

5.1 Boundedness of oscillatory integral operators 

In this subsection we estimate norms of multilinear operators T = JF^" 1 ) defined by (j3.4j) 
and the related composition monomials. The operators JF^" 1 ) have the form ()3.4|) where 
D m = as in (|2~H5]) or D m = [-ir, 7r] d(m_1) as in (jOSJ). The both cases are completely 

similar since we use the same properties of the spaces L\ = L\ {{— ir, n] d j or L\ = L\ (R d ), 

and we do not use in our proofs the boundedness and compactness of the domain [— ir, n] d . 
Hence, we will everywhere consider the periodic case [— vr,7r] d which corresponds to lattice 
equations and without further comment apply the results to the case M. d . 

Lemma 5.1 The operator defined by \3.J$ , \2. 22]) is bounded from E = C ([0, rj , Li) 
into C 1 ([0, r#] , Li) and its norm is estimated as follows 

in 

(u! . ..u m )\\ E < r*Cf" +1 \\x im) \\ J] \\%\\e > ( 5 - 2 ) 

3=1 

\\d T ^ . . .u m )\\ E < Cl m+1 \\ X {m) \\ n W^We ■ ( 5 - 3 ) 

3 
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Proof. According to Condition 12.51 we can diagonalize the matrix exp j— iL (k) ^ j and 
its norm is bounded uniformly in k, T\ and g: 



exp { — iL (k) — 

Q 



< C| Vk g R, g>0,r 1 >0. 



By (HI, JHZD and l|2T22jl 



||^™) (U! . . . y r)|| < Cf" +1 sup x M (k, k) Iff 

|ui (k')| ... iim (k {m) (k, dk' . . . dk^Widk < 

Cl m+1 \\x {m) \\ [ T • • • IK (r!)|| Li dn < r^I m+1 || x (m) || 

./0 



Similarly, 



|a r ^(u 1 ...G m )(-,r)|l < 



Cl m+1 ||x (m) || / / |ui(k')|...|u m (k (m) (k,fc)) 



(5.4) 



• u 



mils 



dk' . . . dk (m_1) dk < 



X 



(m) I 



Ui \\ R .. . U 



■m||£; • 



Corollary 5.2 // M |^JF, T, A, £J (xi . . . x m ) is a decorated composition monomial defined 
by M-ltfy and T is defined by \S. Ji) . \3-J$ then 



M[F,T,\,( ( Xl ...x 



< Cl e+l ri 



niix^iiniw 



E ' 



(5.5) 



1=1 



d T M LF, T, X, C (xi . . . x 



nik (M(7v)) iin 



~ ' I I \\X 

NeT i=i 

where i is the incidence number of the tree T , and e is the number of edges ofT. 



(5.6) 



Proof. We estimate the norm of the monomial M = (Mi . . . M m ) and its time 
derivative applying Lemma I5~T1 Then we use ()5.2j) to estimate II lienor,] L\Y ^ ne f° rma l 
proof is straightforward and uses the induction with respect to the incidence number of a 
monomial. ■ 

Using boundedness of operators we obtain in a standard way uniqueness of solution 
of fl£BD - 

Lemma 5.3 J/ui,U2 G C ([0, To] ,Li) with tq > are two solutions of hS. with the 
same h then ui = U2. 
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Proof. Applying Lemma [4.61 we conclude that 

|| ^ (Ui) - T (U 2 ) \ C {%rr\M) < Cr l \\F - ^ ("2) llcdO.n],^) > < T l < r 0- 

Deriving the above inequality we use that since iV^ < oo the radius in Lemma 14.61 is 
arbitrary large and CV in ()4.13|) according to (|5.2|) is proportional to T\. When the Lipschitz 
constant Gt\ < 1, in a standard way we obtain that ui (r) = u 2 (r) for < t < Ti. Since 
this statement can be applied to Ui (r — T\) and u 2 (r — Ti) we obtain that solutions coincide 
for < r < To- ■ 



5.2 Function- analytic expansion of solutions for modal integral 
evolution equation 

The reduced evolution equation ()3.3|) has the form 

u = .F(u)+x, (5.7) 

where u, x are functions of (k, r). The nonlinear operator T in the right-hand side of ()5.7|) 
is determined by (|3.4|) . x (k, r) = h (k) as in (|3.3|) . We look for the solution operator in 
the form of operator series 

oo 

u = g ( x ) = J2 $ {m) (* (m) ) • ( 5 - 8 ) 

m=l 

The questions related to the existence and the convergence of such series are addressed in 
Theorem 14.71 As a direct corollary of Theorem 14.71 and Lemma 15.31 if applied to the reduced 
evolution equation (|3.3|) we obtain the following theorem. 

Theorem 5.4 Let 

M E <R g = (r.C x <^' +1 )~ 1/(m '~ 1) /8, r* < C^C' 1 . (5.9) 

with C x as in \2.2n}) . C~ as in \2.11 ). Then the series \5. <5j) converges in E = C ([0, r#] , Li). 
The solution operator Q (x) = u determines the solution to \5. 1\) and the operators Q^ m ' 
series &5.8lj) satisfy the recursive relations M.lSj) . 



in 



Proof. iFrom (J2.26j) and (|5.2|) we infer that T defined by ()2.21j) belongs to the class 

A (Cjr, Rjr) if 

r*C x C! m+1 < C T R T n \ m = 2,...,m F . 

If C^Rj} < 1 it is sufficient to verify the above condition at m = mp only. After this we 
apply Theorem 14. 71 where according to (|4.19|) we can take 

R 2 R 2 
Cg= 2{c/ + R?) ' Ra = 4(CrUrY ^ 
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(5.11) 



We take 

C r = Rr= (r,C x C 2 E mp+1 )- 1/{mF ' 1] , C g = 2R g = R^/A 

and apply Theorem [O Note that C S 2 R^ < 1 if r* < C^C' 1 . * 

iFrom Theorem 15 . 41 ( observing that by (|5.11|) R? — > oo when r* — > 0) we obtain Theorems 

Oand Theorem I2~TH1 

To prove Theorem 12 . 1 51 on the superposition principle we apply the solution operator Q to 

a sum of wavepackets h/ (k, 0) as in Definition 12.91 For technical reasons we have to modify 

the wavepackets using cut-off functions described below. 



Cutoff functions. We often use an infinitely smooth cutoff function ^ (77), r\ 6 M. d , satis- 
fying the following relations 

< (77) < 1, (-77) = (77) , (5.12) 
^ (77) = 1 for \r)\ < 7T /2, ^ (77) = for \r)\ > n , 

where 7r < 1 is a sufficiently small number which satisfies the inequality 



< 7r < - mindist {k^, a} . 



(5.13) 



Using \I/ we introduce cutoff functions (k, (3) with support near (k^ defined as follows: 

k- CKi' 



(k,/3)=* 







l-e 



l = l,...,N h . 



(5.14) 



Here e is a small number, 1/2 > e > 0; we take the same e as in Definition |T 
Given a wavepacket h; (k, (5) we introduce a modified wavepacket 

hf (k, 0) = hf + (k, /?) + hf_ (k, /3) , h* (k, /3) = ^ (k, /3) V (k, /3) 

where are defined by ()5.14|) . 



(5.15) 



Proposition 5.5 // (k, (3) is a wavepacket in the sense of Definition \2. iA then hf (k, (3) 
defined by \5.15\) and \5.1J$ is also a wavepacket in the sense of Definition \2.iA and, in 
addition to that, 



h* c (k, / 9) = 0i/|k-Ck*i|>7r o /9 



l-e 



(5.16) 



h, - hf 



< cp. 



Proof. To obtain (l5~TTI) we note that (l2~3"4]) and (HTH^ imply: 



'1 - ^ 



;i-^, c (k,/?))v(k) 



dk < cp, 



(5.17) 



(5.18) 



and ()5.17|) follows. Remaining statements are obtained by a straightforward verification. ■ 
The following lemma shows that we can replace h; by hf in the statement of Theorem 
l2~T5l in particular in ((233), (23BJ). 
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Lemma 5.6 Let satisfy \2. 34 ) ond hf (k, (5) be defined by 15. 15)) . Let 



h/ 



<R, I = 1, . . . , N h where N h R < R c 



(5.19) 



Then the difference 



N h 



N h 



J=i 



Bq,, 



is small, namely 

\\B*\\ E <C{R)(3. 
Proof. Note that since < f ; < 1 we have 



< 



(l-*i,c)hi,< 



< 



(5.20) 



(5.21) 



(5.22) 



and (|5.18|) . Using the Lipschitz continuity of the solution operator Q fsee 14. 6j) and (|5.17|) we 
obtain (15~?T1) . ■ 

Truncation. We will truncate the infinite series (|5.8j) . To this end we define an integer 
m = m (j3 q ) as a solution of the inequality 



2|ln/3 g | < ^ < 2|ln/7 ? | + l 
I In i?c I ~~ lln i?cl 



(5.23) 



where Rg is the same as in ([5.9)1 . We consider then the following partial sum of the expansion 



m=l 

and readily conclude that the following statement holds. 
Lemma 5.7 Let Q be defined by l\5.8\) . then 

Q (h) - Gm(i3) (h) 1 1 < C (R) f3 when h 



(5.24) 



<R<Rc 



(5.25) 



5.2.1 SI-CI splitting for evaluated monomials 

We consider a function h which is a sum of the form ()2.39|) and the solution Q (^F, hj . Ex- 
panding Q( m ^ ( m m ) ) into composition monomials as in ()4.33|) we obtain a sum of composition 



monomials M (J 7 , T) ( h m ) . Then we look at the m-linear monomial M (J 7 , T) (h™ ) where h 
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equals a sum of Nh one-band wavepacket h; as in (|2.39|) . Using the linearity with respect to 
each argument we expand the monomial into a sum of N™ expressions (evaluated monomials) 



M 



(T,T)(^h] = M(r,T)(h h ...h lm )= M (^ T )(u^X ( 5 - 26 ) 

\l=l / h,—,lm h, — ,lm \ i / 



The sum contains evaluated monomials of two kinds: (i) ones which involve the same 
wavepacket; and (ii) one corresponding to the cross-terms (terms involving different wavepack- 
ets). To be precise, we introduce the following definition. 



Definition 5.8 (SI and CI) We say that an evaluated monomial M (J 7 , T) (^h^ . . . hj n 

with the argument multiindexli, . . . , l m G {1, . . . , N} m in the expansion \5.2b]) is self-interacting 
(SI) if 

h = l 2 = . . . = l m . (5.27) 



Otherwise we say that M (JF, T) (h^ . . . h; m ) is cross-interacting ( CI). 
Using this notation we rewrite (J5.26|) : 



N h \ m \ N h 



M(T,T) \\Y>i\ = E M ^> T )((\f) + E M(f,T)(h h ...h lr 

\\l=l J J 1=1 h,...,l m is CI 

(5.28) 

Substituting this expression into (|4.33|) we obtain the expansion 

oo 

g (hx + . . . + h Nh ) = J2 Q ™ + • • • + h ^Y) ( 5 - 29 ) 

m=l 

oo oo 

= E G (( fi i) m ) + • • • + E G ((SD + Gci (> • • • ' ' 

m=l m=l 

where Qq\ contains only CI monomials with cross-terms. 

Proposition 5.9 Every evaluated CI monomial M (JF, T) ^hi, . . . , has a submono- 
mial of the form 

^ [M (J 7 , TO (h h ...h h y..M(F, T s ) (h la . . . h, s ) ) (5.30) 

where all M (J 7 , T x ) ^h^ . . . h^J ; ... ; M (J 7 , T s ) (hi a . . . hi a j are SI, and there are at least two 
indices i and j such that ^ hj . . We call such a monomial a minimal CI monomial. 

Proof. The set of CI submonomials of M (J 7 , T) is finite and it is non-empty since 
M (J 7 , T) itself is a CI monomial. We take CI submonomial of M (J 7 , T) with a minimal 
rank. Its rank is non-zero since every zero rank submonomial is SI. Since the rank is minimal 
all submonomials are SI. Hence it has the form (J5.30|) . ■ 
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5.3 Properties of SI monomials 

According to Definition 15.81 for a SI evaluated monomial we have h/ 1 = . . . = h/ m . Observe 
also that in view of Definition 12.91 every single-band wavepacket h/ has its band number, 
and n' = n" = . . . = nS m \ that is the band rii = uq is the same for all h;. Similarly, 
k*^ = . . . = k„4 m . Having these properties we often omit in this section indices Zj and 
skip n for notational brevity, writing, for example, 

UnX (k) = Wf (k) , u n , c (k) = u c (k) , X ^ = x ™. 
5.3.1 Monomials applied to a single-band wavepacket. 

Here we consider monomials based on oscillatory integral operators and which are applied 
to a single-band wavepacket. We recall that according to (|2.33|) a single-band wavepacket h 
involves two components h + and h and a small complement component hoc. 



Definition 5.10 (frequency matching) We call a decorated composition monomial M I J 7 , T, A, £ 
frequency matched (FM) if for every non-end node N G T the corresponding decorated sub- 

: (m') 



monomial M' = JF^ [Mix' • • • M m , ^{m')J satisfies the following conditions: 

A^oo, C°Voo, j = l,...,m', (5.31) 

and 

m' 

^C W = A, (5.32) 

3=1 



where A, G A defined by &4.4I ), we identify ± with ±1. A decorated composition monomial 



which does not satisfy the above conditions is called not frequency matched (NFM) monomial. 
Collecting separately FM and NFM terms in the expression (j4.51|) we obtain 

M {F, T) (xix 2 . . . x m ) = M T ' ?) ( XlX2 • • • x -) ( 5 - 33 ) 

FM X,C 

+ m(t,tX() (x^.-.xj. 

NFM A,C 

Remark 5.11 Any SI evaluated monomial is either FM or NFM. We do not define for CI 
evaluated monomials if they are FM or NFM. 

Below we show that FM decorated monomials have the following properties which can 
be briefly stated as follows. 
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Property 1. If h (k) is a wavepacket in the sense of Definition 12.91 centered around ik* 
then FM monomial M \^F,T, X,^j ^h m j (k) is also localized about ±k*. This property is 

proved below in Corollary 15.131 

Property 2. The most important property concerning FM-NFM splitting is that the 
result of a NFM monomial application to a wavepacket has magnitude O (g), that is O (j3 ) 
for the scaling ([2.46)1 . Consequently, all NFM terms in ()5.33|) are small (see Lemma 15.161 
below) and they give contribution only to the remainder D in ([2.47)1 . 

Now we formulate exact statements clarifying the above properties. The following two 
statements show, in particular, that an FM monomial transforms a function supported in a 
vicinity of k* into a similar function. 

Lemma 5.12 (operator support ) If Ui,£' ■ ■ .u m ^ m ) are such that 



u m (k (0 ) = when k (0 - C (0 k 



> Si, I = 1, 



and 



k ? = (C + ... + C (m) )k*. 
then JF^" 1 ) ^ . . . u m ^( m )J (k, r) given by satisfies 



k-k 7 



> 5i + . . . + 5 % 



In particular, if the binary indices (, Q m \ are frequency matched (FM), that is 

C = C + ... + C (m \ where ( U \ C = ±1, 
then holds with k^ = £k*. 

Proof. ^From ()3.8[) and ()5.36[) we obtain the equality 

k - Ck* = (k' - C'k*) + . . . + (k< m > - C (m) k,) 

which implies lemma's statement. ■ 



(5.34) 



(5.35) 



(5.36) 



Corollary 5.13 (support of a monomial) If M [J 7 , T, A, £ ) ( hi . . . h m J is a decorated 
composition monomial and 



h L = when 



then 



k (0 _ c w k 
M[F,T,\,C) (h 1 ...h m ) (k) =0 if 



> 8 , I = 1, • • • ,m, 



k- k 7 



> mSo, 



(5.37) 
(5.38) 
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where is defined by \5.3$ . In particular, if M (jF, T, A, Cj (h\ . . . h m j is a FM decorated 



composition monomial, then 

M(V,T,A,C) (hi...h m ) (k) = if |k-Ck*| > mS . (5.39) 

where C satisfies \5.3bS) . In particular, if So = (3 l ~ e and m < Chi (3 then for any 8\ > there 
exists P such that for (3 < (3 we have CV /3 1_e ln/3 < 5\ and 

M (V, T, A, (hi . . . h m ) (k) = when |k - CK I > drop 1 -' In (3. (5.40) 

Proof. To obtain (|5.38|) we apply Lemma 15.121 and use the induction with respect to 
the rank of a monomial. ■ 



Remark 5.14 If M I J 7 , T, A, Cj is NFM and h (k) is a wavepacket localized near ±k*, then 

M (jF, T, A, Cj (h m J (k) is localized near the point k^. As C vary over {—1, l} m such points 
k^ He on a straight line parallel to k*. For m — > oo i/ie closure of the set of such k^ 

a generic k* can fre £/ie entire torus [— ir, ir] d , whereas for the case of C corresponding to an 
FM monomial the closure is just two points ±k*. Hence Property 1 is very useful and, in 
particular, allows to avoid small denominators in coupling terms. 

The following lemma shows that the FM interaction phase function of a single wavepacket 
has a critical point at its center, or, in other words, FM monomials satisfy the group velocity 
matching condition (see [3], [H|). 

Lemma 5.15 If a decorated operator T ™} is FM then the interaction phase function <fi in 

C>C(m) 

hH. <Sj) has a critical point: 

v k0n, c ,n,c (Ck*, t) = att= (C'k*, • • • , C M k*) . (5.41) 

Proof. For FM decorated operator all indices = ± and 

n = n' = ... = n (m) and ( = (' + ... + C (m) . (5.42) 
Hence we obtain from ()3.9|) that 

v k n>CAf (k, k) = cv k c (k) - c (m) v k c (k - k' - ... - k^- 1 )) . 

Since Ck* - (% - ... - C (m-1) ki m_1) = C (m) kl m) and flZHEJ implies 

CV k ^ (Ck,) = C (m) V k ^ (c M Ck*) for C = ±, C M = ±, (5.43) 
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we obtain the desired (|5.41j) . ■ 

Now we consider NFM monomials and prove the Property 2. First we note that (|2.4(J|) 
implies 

u ni (Ki) > > 0, l = l,...,N h . (5.44) 

If = k*, ni — n o satisfy Condition 12. 131 then if ()2.44|) does not hold, ()2.42|) does not hold 
too, hence for m < nip 



m 

i=i 



> lu* > 0, k 



c 



Ec (i) fc* 



(5.45) 



where a;* > is a positive number (we take for notation simplicity the same small enough 
const sant in ()5.44|) and f)5.45D . 

The following Lemma, which is a version of the standard statement of the stationary 
phase method, shows that the action of an NFM monomial on a wavepacket produces a wave 
of a small amplitude. 

Lemma 5.16 Let the decoration projections be defined by \5.1)) . Assume that Condition \2.13\ 
holds. Let indices . . . , Q be NFM, that is either one of them is oo or 



C^C' + --- + C M , c W) = ±i, C = ±i- 

Let S nfm > be small enough to satisfy 



(5.46) 



5 NFM max \Vui (k)| < -u+, 1 = 1,..., N h , 

|k*i-k| <S NFM 4 



where cu* is given in fl<5.^<5| ). Let k,k^ satisfy \3.12\) and be such that 



(5.47) 



E|k^-c (j) k* 

3=1 



< s 



NFM, 



k-k. 



< s 



NFM, 



(5.46 



where k^ zs defined by \5. 34\} and k* = k*i satisfy the conditions \5.44\ ) an d \5.4ty - Let the 
functions u- ^ (k, r) satisfy the condition 

Uj^j) (k, r) = when = oo and u.^tj) ^C^k* + s, r j =0 when |s| > S^fm- (5.49) 
T/ien 



C,C',-,C (m) V Ul ' C ' ' ' ' U ™.C (m) 
2^r 



<l£|| x M|| C 2m+iTTi| G 
/.i.ii " J. J. 



\E UJ 



■j \\e 



(5.50) 



^ II II ^J.J. 
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Proof. If one of the indices C M equals oo by (lo~473|) (m) = and (l5~o7H) is 

satisfied. Now we consider the case when all (, . . . , Q m ' are finite. We denote for brevity 
Lu no = lu, k*i = k* and n ^ - ^ = <fi. Since (|5.48jl holds we get from ()3.9j) that 



> (k, fc) - (k, | < \u (k') - w (c'k,; 

m 

< max I Vcj (k) I V \k U) - ( ij) k 4 
|k*-k|<,w/ 



a; 



(k^) -u (C (m) k, 



Using (|5.47|) we conclude that 



k,Jfe 



> 



k, k* 



< Snfm max |Va>(k)| 

|k» — k\<5rtFM 



4 1 *l 



(5.51) 



By (|5.4fjjl the condition ()2.44j) is not satisfied, therefore (|5.45jl holds and implies that 

>lu«. (5.52) 
Using (l532"j) . flOg) and (l5~4"TI) we conclude that 



k^, 



k, k* 



> lu* 



lu (k) — lu I k? I > c<j* — 5atfa/ max | Vlu (k) | > -a;*. (5.53) 

v " ; lk*-k|<(5jvFM 1 



4 



Together with (j5.51|) this inequality implies that when (|5.48J) holds we have the estimate 



k,k 



1 

> —U)t 

~ 2 



(5.54) 



Now we note that the oscillatory factor in ()3.8|) 



exp \ \(f) ( k, k ) — 



5 J i0 (k, fc 
Integrating (|3.8jl by parts with respect to T\ we obtain 

7 (m) 



<^ r , exp <j in f k. /»' ) — 



f?exp > 


\<j> (k,kj 


5} 




i> (k, fc) 





^7(u 1 ...u m )(k,r)= (5.55) 
(k, fc) fi 1>{ , (k\ r) . . . £ m>C ' (kH (k, fc) , r) d^ d £ 



£ (m) 



id (k.k 



X ( ™> (k, fc) fi 1>c , (k\ 0) . . . u mX , (k™ (k, k) , 0) d^ 1 )^ 



i(/> [ k, k 



exp < i(f) ( k, 



Q 



X 



(m) 



m) 



k,k)d T1 \u ljC {k')...u mtC (is.™ (k,k))] ~S m ~ l Hd Tl . 
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Estimating the denominator by (|5.54j) and using (|3.5j) we obtain (|5.5U|) . Finally, we consider 
the case when £ = oo and all remaining indices equal ±. We expand into sum of 
IL^ as in ()4.38|) . In this case \^ (k, involves a projection H n ,( with n ^ «o (the 
oscillatory integral may involve — 1 terms with such n). For a fixed n the corresponding 
phase function ^k, k^j takes the form 

(k, k) = <f> n ^ (k, k) = (u n (k) - ('u no (k') - ... - C^u no (k^) . 

Using again ()5.45|) (now with n ^ no) we obtain that ()5.52|) holds. This implies ()5.54|) as 
above provided Snfm is small enough. Hence, the relation ()5.55|) holds, implying readily the 
desired bound ()5.50|) . ■ 



5.3.2 FM and NFM monomials for SI oscillatory integrals 



The following below theorem shows that NFM monomials are of the order O (g) as g — > 0. 
We begin first with the following statement. 

Lemma 5.17 Assume that Condition \2. Up holds. Let a monomial S = J 7 ^ . . . M s ^(s) 

have all submonomials M l ^ . . . M s which satisfy FM condition \5. 3b]) . but S itself is not 
FM. Assume that S is applied to wavepackets h/ which satisfy Definition \2. M and 



(k,/3) = 0z/|k-Ck^|>7r /? 1 - e ,C = ±. 



(5.56) 



Then 



< 



4p \\x 



Ml 



\uj (k* 



c i s+i nl! M /. 



+ 



(5.57) 









w(k»)| 



c E +i j2\\ d ^M\ nlKc^ , e=c([o,tj,l 1 ). 



i=l 



Proof. Since M l ^ . . . M s are decorated FM submonomials we can use Lemma [5.121 
and Corollary 15.131 Applying Corollary 15.131 and using ()5.12|) we obtain that 

M ic(i) (k w ,n) = when k (0 - C (0 k* > u (M hc( iA f3 l ~ € -K Q ,l = 1, . . . , s. (5.58) 

where v (M) is homogeneity index of M. Consider now the oscillatory integral (J3.8j) which 
determines S, namely 



^(M liC(1) ...M S)C(s ))(k,r 



exp < l 



k,k)- 



(5.59) 



X 



» 



k, k) M 1>f( i) (k', n) . . . M S)C(S) (k« (k, fc) , n) d^kdr,. 
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We apply Lemma 15.161 where, according to (|5.58|) and (|5.62|) Snfm = "mp e vr . According 
to (I53TI1) 



\\si 



+ 



C,C(a) 

4^r, |M s) 



M 1jC (d . . . M S)C ( S) ) (k. 7-; 



(5.60) 



that implies 1)5. 57j) . 



-Cf +1 \\drM iiCW II J| ||M. if0 ) ,E = C([0,t.], In) 



i=l 



Theorem 5.18 Suppose that (i) the inequalities \5.44 ) hold; (ii) h; are wavepackets in the 
sense of Definition \2. 9i (Hi) the relations \5.5b}) hold; (iv) the projections are defined by \5.1\) ; 
(v) Condition ^. lc\ holds. Then a NFM decorated monomial based on oscillatory integrals T 
defined by \3.4D satisfies the estimate 



Aqt\ 1 [l-ffl J r 



M{F,T,\,() (h!...h 

re Yl 



< 



X 



c([o,T,],ro 

m 

""■^'11 J]||h 



(5.6i; 



NeT,r{N)>0 



1=1 



where i, m and e are respectively the incidence number, the homogeneity index and the number 
of edges ofT. 

Proof. Let M (^F, T, X^, C( m )j ■ ■ • h m j be a NFM decorated m-homogenious mono- 



mial. We find a decorated submonomial S = M {J 7 , T (N ) , A( g ), C(m)J °f M yJ 7 , T, A( 9 ), C( m ) 
with such A^o that S is NFM and has minimal rank of all NFM submonomials. We denote 
by ro the rank of S, by i' its incidence number and by s = v (S) = v (T (No)) its homo- 
geneity index. This monomial has the form S = J 7 ^ (m x ^v, . . . M s V Since the rank is 

minimal, all decorated submonomials M, m) are FM and their ranks do not exceed ro — 1. 
Then according to (J4.21)) their homogeneity indices satisfy 



v (m ic(1) ) +... + V {M s& ^j = s < 



m. 



(5.62) 



Applying Lemma 15. 171 we obtain (|5.57|) . Now we use Lemma I5~T1 and 15.21 Applying Lemma 
we obtain 



\sh n 

NET\T'(N ),r(N)>0 



M({^},r,r)(h 1 ...h. 



< 



tt(MJV)) tt r 



n 



l>x+v(T'(N )) 
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Note that the norm of 



T(N) 



is estimated by ()5.2|) and norm of S by f)5.57j) . In turn, we 

estimate right-hand side of ()5.57|) using ()5.2|) and (|5.3|) Taking into account that s < m in 
the sum in (|5.60j) we get the estimate (|5.6 lj) . ■ 

We also consider the case when Condition 12.131 does not hold and Condition 12.231 holds. 
In this case we give an alternative definition of FM and NFM decorated monomials. 

Definition 5.19 (Alternative Frequency Matching) We call a decorated composition 
monomial M (^F,T,X,(^ alternatively frequency matched (AFM) if (i) every node of T 
has an odd number of child nodes (at least three); (ii) for every non-end node N G T the cor- 
responding decorated submonomial M' (^F, T (N) ,X,Cj = ^ i^M\ £ ■ ■ ■ M m , ^( m ') j satisfies 
/ TP71) and 

sign(^2C Q) j =A, (5.63) 

where A, G A defined by d^.^i| ), we identify ± with ±1. A decorated composition monomial 
which is not AFM is called alternatively not frequency matched (ANFM) monomial. 

Now we prove a statement analogous to Theorem 15. 181 when Condition 12.231 holds. 

Theorem 5.20 Assume that assumptions of Theorem NFM hold with Condition \2.1,Si re- 
placed by Condition \2.2'A Then \5. holds. 

Proof. According to Corollary I5.13[ if h; 1 = . . . = h; m = h; satisfy Definition 12.91 and 
then M (>, T, A, () (h h . . . h im ) = M (>, T, A, f) (h,^ . . . h Jmif ( m) ) has support in 
a m(3 l ~ e vicinity of = z/k* defined by ()5.34j) . v and m are odd integers, m > 3, v — 
C + ... + C (m) - Let S = M (V, T', A, f) be minimal ANFM submonomial of M (V, T, A, , 

that is if T" C T' then M (^F, T", A, (fj is AFM submonomial of S. The monomial 5* has 
the form of (J5.59)) with the interaction phase function 



(k, k) = (uj n (k) - C'^no (k') - ... - C (S W, (k (s) ) . (5.64) 



The integrand is non-zero near k^ = i/jk^and applying ()5.63j) to every AFM submonomial 

we get 

C (,) = sign (i/,) . (5.65) 

Using (EUSD and (ITTH|l we obtain 

4> c ^ (^K, K^j = (u n (vK) - C'u no (z^k*) - ... - ( {s) u no (v s K) (5.66) 

= (u n (vK) - sign (u{) (k») - ... - sign (u s ) \u s \uj no (k*) 

= C \v\ uj n (k*) - (z/i + . . . + v s ) u no (k*) , v = v x + . . . + v s . 



62 



Therefore, since S is ANFM, ( ^sign(z/) and since v is odd, 



^z/k*, j = —2vuj na (k*) 7^ 0, (5.67) 



therefore (|5.52|) holds. We can repeat the proofs of Lemma 15 . 1 61 and Lemma 15. 171 and obtain 
1)5.57]) . From 1)5.57)1 we obtain 1)5.61)1 as in the proof of Theorem 15. 181 ■ 

Below we give estimations for the derivatives with respect to k of a composition monomial 
applied to a wavepacket. Note that ([2.35)1 admits a singular dependence on (3 of wavepackets 
h^(/3,k). This type of dependence also naturally comes from explicit formulas as ()2.36j) 
which yield that the first derivative with respect to k has a factor @~ . Below we estimate 
dependence on (3 of monomials applied to wavepackets and will show that they have the same 
type of singularity. 

Observe that by ()5.13|) all the points k^ are at the distance at least 27Tq from a. Hence, 
according to Definition 12.31 and 1)2.28)1 

max (|V 2 k u,| + |V k u,|) < C u * (5.68) 

|k±k t (|<7r , l=l,...,N h , 



max 

|k±k» ; |<7T , l=l,...,N h 

The following seemingly technical Lemma describes a very important property of solutions. 
It shows that the k-gradient of solutions behaves, roughly speaking, as the gradient of initial 
data. Corresponding estimates play a crucial role in the control of smallness of interaction 
of different wavepackets. 

Lemma 5.21 Let M (^F, T, A, (f\ ^h^ . . . hi m ~\ be a decorated monomial which is SI. Assume 

that hij = hi are wavepackets satisfying Definition \2. yi \5.5b)) and \5.1ty) . that \2.4b\j holds 
and 

^-"m < tt . (5.70) 

Assume that either Condition \2.1cH holds and the monomial is FM or Condition \2.2<A holds 
and the monomial is AFM. Then 

|v k M (>,T,A,C) (h h ...h lm ) 

where E — C ([0, r*] , L-y), r* < 1, with i = i (T) and e = e (T) being respectively the incidence 
number and the number of edges ofT. 

Proof. We use the induction with respect to the incidence number % of a tree T. First, 
we consider the case when Condition 12. 131 holds and M (^F, T, A, Qj is FM. For i = ()5.71j) 
follows from ()2.35)) . Now we assume that ()5.71]) holds for the incidence number less than i 



V X ^(k,k',...,k( m )) 



< c x ci 



m+1 



(5.69) 



< CC x T\Ci +e C l - l R m - l l3- l - e m 2 , (5.71) 
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and prove it when the incidence number equals i. Since arguments of M (^F, T, A, (fj are SI, 
according Definition 15.81 U — . . . — l m — I. It is sufficient to prove the boundedness of 



V k MLF,T,A,C [hi 



V k J^ s) [M UI ...M S 



where .\/ : . . . M s are decorated submonomials, Mj^ = Il^Mj. Let the submonomials have 
incidence numbers i\, . . . , i s and homogeneities m 1; . . . , m s respectively satisfying 



%x + . . . + i s = i — 1, mi + . . . + m s — m. 



We have by fl£B) 

V k ^ s) (M 1)f ,...M a>cW )(k,r) = V k 




J[—ir,w 



(s-l)d 



exp < l 



By Leibnitz formula 



where 



(k, fc) M u , (k') . . . M SjC ( s) (V s > (k, fc)) d^-^^jfcdri. 
V k ^l s) (m 1i? , . . . M, iCW ) (k, r) = A + J 2 + J 3 , 




7"! 



V k exp < icj) , ? k, k 
^[-tt.tt]^- 1 ^ l I ,4 v 7 e 

(k, fc) m 1jC , (k) . . . M S;CW (k« (k, £)) d^-^^dn, 



(5.72) 



(5.73) 



(5.74) 




'- 2 : / / exp < 10 A * ( k, k ) — 

v kX H (k, £)] m 1jC , (k) . . . M s>&) (k« (k, £)) d'-^foZn, 




*/ [ — 7T,7T 



(s-l)d 



exp < l 



k,fc 







to 



k, fc) M lif , (k') . . . V k M sC ( s) (V s ) (k, d (s - 1)d A?rfr!. 



By dH3D 



M jxU) (k&'>) 



< CT i(3) cf J)+eU) cfR m ^] = 1, 



(5.75) 



Using (|3.5|) . (|5.75j) . (|5.72|) and the induction assumption we get 



s-1 



i/ 3 i< iix (s) iinii M ^ ( kW ) 



V k M SjC ( s) 
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t*ri < CC?K n - 1 T i C!g +e C i x l3- 1 - e . (5.76) 



From ()5.75|) and the smoothness of x ( k, k ) we get 



|/ 2 | < c , /3- 1 -v i c 1 m c| i+e c; J R m . 



Now we estimate Ji. Using ()3.9|) we obtain 



(5.77) 



h 



eX P 1 ^ A,C I k ' fc 



Q 



(5.78) 



A.C 



AV k cj (k) + C {s) V k u (k (s) (k, k^ 
k, Pj Af 1>c / (k') . . . M SjC ( s) (k (s) (k, fc)) d (s - 1)d ^dri. 



The difficulty in the estimation of the integral I\ comes from the factor ^ since g is small. 

Note that according to p 2 /g<C. Since M ^JF , T, A, is FM, its every submonomial 

is FM too and we can apply to them Corollary I5.131 which yields 



M u(j ) (k^) = for - C U) K 
Hence, it is sufficient to estimate iifor. 



> mjiioP 1 \ j = 1, . . . ,s. 



k (i) _ <Ki) k< 



< 5\ = nm f3 6 for all j. 



V k A ^ ( Ak*, k* 



0. 



(5.79) 



(5.80) 



According to Lemma [5.1 51 since A, ( are FM 

-AV k cu (k*) + C (s) V k cu ((k« (k„ k 
Using ()5.68|) we conclude that in a vicinity of k* defined by (|5.79j) we have 
-AV k cu (k) + C (s) Vk^ (k (s) (k, jfe))l I < 2 (s + 1) C^. 
This yields the estimate 

\h\ < CCi i+e T i C i x C^ l (3- 1 - e m 2 R m . (5.81) 

Combining (|5.8ip . (|5.77|) and (|5.76j) we obtain ()5.7H) and the induction is completed. Now we 
consider the case when Condition 12 .231 holds and the monomial is AFM. Note that according 
to Corollary 15.131 the submonomials M-^) have supports near z^k*, with an odd Uj. By 
Lemma [5.121 the monomial itself is non-zero near z/k*, v = v\ + . . . + u s ; since s is odd v is 
odd too. Obviously, one of Uj has the same sign as u, we assume that j = s, that is 



sign [y s 



sign [Pi + . . . + v s 



sign \y) 



(5.82) 
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the general case can be reduced to this by a relabeling of variables. The interaction phase 
function is given by (|5.64|) and since the submonomials are AFM (|5.65|) holds. According to 
(BUB} V k (w (-k)) = - (V k w) (k). Therefore, using (ETTIoT) we obtain 

V k A ^ (Vk*, k*j = AV k w (i/k*) - C (s) V k ^ (i/ s k*) = 
A (Vkw) (sign (i/) k*) - C (s) V k u; (sign \u s \ k*) 
= A (V k w) (sign {y) K) - C (s) V k w (sign (v 8 ) k») = (Asign (i/) - C (s) sign (z/ s )) (V k w) (k*) . 

Using (|5.65j) we conclude that 

V k A ^ ^k*, fc*^ =0, k* = (i/ik*, . . . , z^sk*) . (5.83) 

Using ()5.83|) instead of ()5.80|) we conclude as in the first half of the proof that (j5.71|) holds 
in the AFM case too. ■ 



5.4 Properties of minimal CI monomials 

Here we consider CI evaluated monomials with arguments involving different wavepackets h;. 
Since the group velocities of wavepackets are different, namely (j2.41j) is satisfied, there exists 
Po > such that 

|Vc^ (K h ) - Vcu (k, /2 )| > po > if l x ± l 2 . (5.84) 

The next lemma is a standard implication of the Stationary Phase Method in the case when 
the phase function has no critical points in the domain of integration, namely when ()2.41|) 
holds. 



Lemma 5.22 Lefk^i anduj n be generic in the sense of Definition ^. 2\ Let T 1 ^ be defined 
by \3.4]) , rn (0) be as in \5.23\) . We assume that \2.2ffy and \2.41]) hold. We also assume that 
Urm . UTZ(\) . JPfl J, fOSp and WJ®> hold. We assume that M (F,T) (h h . . . h lm ) is a 

monomial with homogeneity index m evaluated at arguments with CI multiindex Z 1; . . . , l m , 
but every evaluated submonomial of M (J 7 , T) (h^ . . . hi m J is SI. Then for m < m{f3) and 
small (3 



M(F,T) (h h ...h 



< T 

E Po 



m 2 R 



m—l 



(5.85) 



where i and e are respectively the incidence number and number of edges of T , R is as in 

Proof. Since k^ are not band-crossing points, the relations (|5.69|) and ()5.68|) hold. We 
expand M{T,T) into a sum of decorated monomials M [J 7 , T, A, () as in (j4.51|) . which 
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contains no more than 3 l ( T )+ m terms, and i (T) + m < 2m. The submonomials of every 
decorated monomial are SI by the assumption of the theorem. If Condition 12.131 holds, the 
submonomials are either FM or NFM; if Condition 12.231 holds, the submonomials are either 
AFM or ANFM. If a decorated submonomial M LF, T', A , ( ) is NFM we use Theorem ETTBI 
and obtain from (J5.61|) the inequality 



m (jf, r, a , c ) (h v+1 • • • K^)\\ B ^ t 1 + m l ci'^'c^ 1 



(5.86) 



where %' and e' are the incidence number and number of edges of the subtree T' . Alterna- 
tively, if Condition 12. 231 holds, and .a decorated monomial M (^F, T', A , C j is ANFM, we use 

Theorem I5.2UI and obtain from (J5.61)) the inequality (J5.86j) . Using (|5.86|) in both cases we 
obtain 



M F,T,A,C [h h ...h 



< CQrt 1 Ci +e C i x mR m . 



(5.87) 



Now we consider the case when Condition ^. 13l holds and every submonomial of M yj 7 , T, A, ( 

is FM. We write the integral with respect to T\ in (|5.59|) as a sum of two integrals from to 
(3 and from (3 to r, namely 



^ s l(M 1 ...M s )(k,r)=F 1 + F 2 , 



(5.88) 



exp < l 



k, jtj ^ | A^l (k, fc) K-Wkdn, F 2 = £ ... dr x 



where 



A (s l ( k, k 

CiC(m) 



x;, lk,£)M 1 (k')...M s (kM) 



(5.89) 



Mj are submonomials of M. According to Corollary 15.21 with r* = /3 

II^IL, < 2C| +2s C x /3n \\M 3 \\ E < *V 2 V f J[ ||h,. < pC^t^R™. (5.90) 

1 1 

i=i i=i 



Now we estimate F x . Since M(jF, T) is CI, there are two SI submonomials M~ 1 and M, 



applied to (h^ j and (h^ 



mi 



1112 



.12 



with 7^ lj 2 . Let us assume that lj 1 = h, lj 2 = l s (the 



general case can be easily reduced to it by a relabeling of variables). We denote 

<t>' = V k /0 Ci? (k, t) = Vuu (K h ) - V k(s) ^ (Ki a ) ^ 0, p = 07 \<j)'\ . (5.91) 
By (jEHlD and (jEHD we obtain 



p-V k '0 Ci ^ ( k, t 



> po > for k = k* = (k* tl , 



k *0 



(5.92) 
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Note that 



exp { i0 ( g ( k, k ) — - 



Q 



Q 



-p V k ' exp <H0 ? k, fc — . 
ip-V^ a [k,k)ri I ^ V J B 



Tl 



Using this identity, ()2.25|) and integrating by parts the integral which defines F\ in ()5.88|) 
we obtain 

Fi = J T I (k, n) dr 1} I (k, n) = jf exp |i0 c ? (k, kj ^ j (k, fc) d*- 1 ^ = (5.93) 



g exp (i^ ? f k, fc) ^ ) AH ( k, k ' 



V k '0 f ^ ( k, fc ) ■ p 



From (|5.56jl . Lemma 15. 121 and Corollary 15.131 we see that in the integral I (k, t-\) the inte- 
grands are nonzero only if 



< TrijTioP 1 e , |k — ^k*| < mTToP 1 e ,m 1 + . . . + m s < m, 



(5.94) 



where ir < 1. Using the Taylor remainder estimate for 0^ at A;* we obtain the inequality 

V k '0 Cj £ (k, fc) - 0' < Sm/? 1 "^^^ if (|5IMD holds. (5.95) 
Suppose that (3 is small and satisfies 



Po 



(5.96) 



Condition ()5.96|) is satisfied for small /3 if m < m(/3) as in ()5.23|) . Using ()5.95|) we derive 
from (lo~m (IBUBl and lETSEI) that 



p-V k '0 c ^ (k, /c 



> y > if dniHD holds. 



(5.97) 



Now we use (|5.97|) to estimate denominators, (|5.68|) to estimate second k'-derivatives of u 
and ()5.69|) to estimate Vk'X- We conclude that 



|/(k,T!)| <Cf +] 
Q 



TlPo 



V k ^g(k./, 



+ 



8gQ 



£J,2 



A$ (k.A- 

t.;t,(m) 



d (s - 1)d fc (5.98) 



< 



(V*-V fcW )x« (k,-)|| + ^||x (m) (k 



8a. 



TlP0 

(k, 



2 IM m ) 



||M,|| Li \\VvM 1 \\ Li + ||M,|| Li ||V k(s) M s || Li 

J=2 J=l 
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To estimate VMj we use Lemma 15.211 We also use (|5.2|) and (|5.5|) to estimate HM,-^ . 
Therefore, using (|5.72|) . we obtain 

|/(k,n)| < ^r^Cl^C^-^m 2 ^- 1 . (5.99) 

Finally, we consider the case when the alternative Condition E2S1 holds. In this case Mi and 
M s according to Lemma 15. 121 are localized near fik*^ and V2^*i a with some v\ and v-i, we 
use fl2~Ml) to obtain both for AFM and ANFM submonomials. Therefore flBHTfl holds 

and we again get ()5.98|) and ()5.99|) . So, we proved ()5.99|) in all cases. Integrating ()5.99|) in 
Tx we obtain 

\\F X \\ E < Cr\- l C 2 ^C l x -^m 2 |ln/3| R m ~\ (5.100) 

P Po 

Using summation over all A, ( (the sum involves no more than 3 2m terms) we obtain ()5.85|) 
from flUEZD and fTKKty . m 

6 Proof of the superposition theorems 

In this section we prove Theorems 12.151 12.191 on the approximate modal superposition prin- 
ciple. 



6.1 Proof of the Superposition principle for lattice equations 

Here we prove Theorem 12.151 First we note that according to Lemma 15.61 we can replace 
hj by h* in the statement of Theorem 12.151 in particular in (|2.47J) . ()2.48|) . Hence we can 
assume that ()5.56|) holds. 

Based on Theorem 15.41 we expand the solution of ()2.Hj) into series (15.8)1 and then into 
the sum of composition monomials M (J 7 , T) as in (j4.Hrij) : 

oo 

g (>, h) = h + 53 °? M T ) (h m ) . 

m=2 TgT m 

where 

N h 

i=i 

and the relation (|5.19j) (that is N h R < Rg) holds, where Rg is the radius of convergence 
from Theorem 15 A\ R will be specified below. Using Lemma 15.71 we conclude that 

m(/3) 

S(V,h)=h+£ C ? M T ) (h m ) +9,\\9\\e<P, ( 6 - 3 ) 

m=2 TeT m 



< R, I 



N h , 



(6.2) 
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where m (/3) is denned by (|5.23p . Then we expand every monomial M (JF, T) yh m j according 

to ((ElHI) into the sum of the terms M (J 7 , T) (h h . . . h z A Note that since m (/?) < C \ln(3\, 
conditions (|5.96j) . (|5.7U|) . (|5.47|) are satisfied if (3 is small enough for every m <m {(3). The 
monomials M (J 7 , T) (h.^ . . . h^ m j belong to two classes, SI and CI (according to Definition 
15. 8|) and the class is determined by the multiindex (Z 1; . . . , l m ) = I. Using ()6.3|) we conclude 



that 



V 1=1 J 1=1 

6 = EE E c T M(^T)(h, 1 ...h u )+^ |M| B <C/3. 



(6.4) 



m(0) 



m=2 T£T m CI ii,...,i„ 



To obtain ()2.48|) we have to estimate the sum in D and show that it is small. It follows from 
flP3)) that 



m(/3) 



EE E ctm^t)^...^ 

M(F,T) (h h ...h t 



m=2 T<=T m Clh,...,l 
m(f3) 

< E ^ E cr su p 



< 



m=2 TeT m reT m ,GH 
m(j9) 

<^iY™ CoC ™ sup 



< 



m=2 



TeT m ,CT 



M(.F,T) h ;i ...h 



Now we consider an evaluated monomial M (JF, T) (h^ . . . h/ 



with arguments given by CI 

multiindex / = . . . ,l m ). To prove that this monomial has a small norm, according to 
Lemma 14.241 it is sufficient to show that one of its submonomials is small and the relevant 
operators are bounded. According to Proposition 15.91 the monomial M (JF, T) ( h; 1 . . . hi 



contains a submonomial M (JF, T') ^h; a , . . . h^„ J with the homogeneity index s = s" — s' + 1, 
the incidence number i' and the rank r' which is minimal in the following sense. The monomial 
M (JF, T') (h-i s , . . . hi s „ j is CI, but every its submonomial M (JF, T") ^h^„ . . . h^,,, j is SI. Now 
we use the space decomposition (|5.1|) and expand M (JF, T") as in (J4.44j) into a sum of no more 
than 3 2m decorated monomials M (^F, T', A, ( j (h^, . . . h^„ j . The decorated submonomials 
of every decorated monomial are SI. We apply Lemma [5.221 and conclude that 



m(f,T'X<) (hL, - h 



< C 







l+e 



|ln/3|+/3 



-ri- i d +2if 

Po 



CIR S 

A, 



(6.5) 
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Hence, there is a submonomial of M(J-,T) (h^ . . . h; m J with a small norm. Namely, since 

()2.46|) and ()2.5|) are assumed, this small submonomial provides the smallness of the norm of 

the whole monomial M (J 7 , T) (h.^ . . . h; m ^ according to Lemma \A. 241 We also use Corollary 

15.21 and ()2.26|) to estimate norms of remaining submonomials of rank r and apply ()4.32|) and 
(f5~72"J) to obtain 



M{F,T) (h h ...h 



< 3 a 



2 m 



P 



1+6 



\hx/3\+p 



Cim 2 T i - 1 C^ +2i C! ( R m - 1 . (6.6) 



Since e = % + m — 1, using ()4.25|) and the inequalities z (T) = i > m/mp, i < m — 1 we get 



m(/3) 



EE E c T M(F,T) (h h ...h h 



m=2 TeT m Clh,...,l r , 



(6.7) 



< Co 



fc|ln/?|+/? 



J m=2 



with ci = 9C|C X . The series converges if, in addition to ([5. 19)1 . R satisfies the inequality 

RN hCl T l J m r < 1. 

For such R and r#, combining (|6.7|) with ()6.3|) and using ()2.46|) we obtain (|2.48|) . and the 
Theorem 12. 151 is proved. 



6.2 Proof of the Superposition principle for PDE 

Here we prove Theorem l2.25l fand its particular case Theorem 12. 19)) . The proof is completely 
similar to the above proof of Theorem 12.151 up to every detail. One only have to replace D m 
given by ()2.23|) by D m given by ()2.65|) and the space L\ is now defined by ()2.66|) instead of 

(USD- 

Remark 6.1 Note that smallness of CI terms is essential and is based on different group 
velocities of single band wavepackets. Note that separation of different wavepackets based 
only on FM and NFM arguments as in Lemma 5.11 is impossible since there are always FM 
monomials with different I because of the symmetry conditions \2.15\) . \2.1b\) . for example 
FM condition 



C^n,C (C k *) - CW (C' k *l) - C^n" (C" k *2) - C"W" (C'" k *3) = 



is fulfilled if 



n = n\ C = C, K = Ki, n" = n"\ (" = -("', k, 2 



>-*3 



independently of the values of k* ; k*3 and independently of a particular form of functions 
LU n (k) . 
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7 Examples and possible generalizations 

7.1 Fermi-Pasta-Ulam equation 

FPU equation on the infinite lattice has the form 

d 2 t x n = (x n -i - 2x n + x n+ i) + a 3 ((x n+ i - x n f - (x n - x n -i) 3 ) 

+«2 {{x n+ i — X n ) — [x n — X n _i) ) . 

It can be reduced to the following first-order equation 

9tX n — Un ~ Vn-li &tVn = %n+l ~ %n + ^3 (x n+ \ — X n ) + OC2 {x n +l — 

We introduce lattice Fourier transforms x, (k) and y (k) by (|2.2jl . namely 

(k) =J2^e-' mk , kE [-7T, 



(7-1) 



(7.2) 



,7T . 



First we write Fourier transform of the linear part of (|7.2|) (that is with 03 = a 2 = 0). 
Multiplying by e~ and doing summation we obtain 

$x (Jfe) = y (k) - e- ifc y (A;) , d t y (Jfe) = e ik x (k) — x (k) . 
that can be recast in the matrix form as follows 





X 




. V . 





e ik - 1 



(e ifc -l) s 






X 




. y . 



The eigenvalues of the matrix are purely imaginary and equal iu^ (k) with 

k 



U((k) = C e 



iA- 



2( 



sm 



( = ±, — 7T < k < 7T. 



The eigenvectors are orthogonal and are given explicitly by 

sc(*) 



1 


" i( \e ik - 1 " 


1 


< 

e ifc -l 


y/2 \e ik - 1| 




~V2 



C = ±,k^o. 



(7.3) 



Now let us consider nonlinear terms. Note that the lattice Fourier transform of the product 
x (n) z (n), n G Z d is given by the following convolution formula 



xz ( k ) = 7^~\d I x(s)z(k- s) ds 

27T J[-7r,7rl d 



(2 

as in the case of the continuous Fourier transform. Note that 

x n+1 - x n (k) = (e lk - l) x (k) , 
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(7.4) 



and, hence, the Fourier transform of the cubic term of the nonlinearity in (|7.2j) is 



(7.5) 



(2tt) 2 Jk'+k"+k'^=k 



Ak' 



1 ) ( e ifc " - 1 ) ( e lfc "' -1)5 (k') x (k") x (k'") dk'dk 



Ak" 



and similar convolution for the quadratic term. 



7.2 Examples of wavepacket data 

Here we give examples of initial data for PDE in IR d and on the lattice Z d which are wavepack- 
ets in the sense of Definition 12.91 We define a wavepacket by ()2.33|) where is chosen to 
satisfy (|2~85|) and 

Recall that a Schwartz function is an infinitely smooth function $ (r), r 6 M. d which 
rapidly decays and satisfies for every s > the inequality 



sup |rr|a?$(r)|t2r<Ci(s) 



(7.6) 



where 



|«|+p<s 



9 r Q $ (r) = d" 1 . . . $ (r) , a = 



, a d ) , a = oti + . . . + ad- 



It is well known that Fourier transform of a Schwartz function remains to be a Schwartz 
function and that its derivatives satisfy the inequality 



sup l k l P ^( k ) <C 2 (s). 

k |a|+p<s 



(7.7) 



Example 1. We consider equation in R d as in Subsection 1.2. The simplest example of a 
wavepacket in the sense of Definition 12.91 is a function of the form 1)2.36)1 where 



h c (k) + V k Mk) 



+ Ik 



h c (k) 



dk < oo, . 



(7i 



and g n ^ (k) is an eigenvector from ()2.13)) . Note that (3 d h^ (k//3) is the Fourier transform of 
a function (/3r). 

Lemma 7.1 Lei h k) , k <G M d 6e defined by \2.3b\) . \7.$ . Then k) is a wavepacket 

with wavepacket center k* in the sense o f De finition \2. .91 with L\ = L\ (W 1 ). 



Proof. First, condition ()2.32j) holds since 

'k-Ck* 



h c (ft 







-d 



he 



5n,C ( k *) 



dk 



(k*) 



dk. 
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Condition (|2.33|) is obviously fulfilled since 

h c (/3,k) = n n , c (k) h c (ftk) 

Inequality ()2.34j) follows from the estimate 



P' 



\k-(k,\>p L 



ht 



k-Ck. 

P 



dk<p I |k 

|k|>/3" e 



fcc(k) 



dk < C/3. 



(7.9) 



To verify (|2.35j) we note that since H n ,c (k) smoothly depend on k near £k* we have 



/ 



lk-Ck.l^/3 1 



<<7 / /? 

f |k— Ck, |</3 1 ~ E 



-d-l 



V k /i, 



V k h c (/3,k) 
k-Ck* 



dk 



P 



< cp- 1 [ V k /i c (k) 



+ p- d 

dk + C 



hi 



k-CK 

P 



dk 



and (JZIBD implies (l2~351) . 



Example 2. Let us consider a lattice equation in Z d as in Section 1.1. We would like to 
give a sufficient condition for functions defined on the lattice which ensures that their Fourier 
transforms satisfy all requirements of Definition 12.91 We pick a Schwartz function $ (r) (see 
(I7.fjjl ). a vector k* e [— 7T,7r] d and introduce 



h(p,r) = e" irk *$ (J3r) , r e M d . 



(7.10) 



Then we restrict the above function to the lattice Z d by setting r = m. The following lemma 
is similar to Lemma f7. 11 

Lemma 7.2 Let $ (r) be a Schwartz function, h^(P,r) be defined by \7.1(A) , h^(P,k) 
be its lattice Fourier transform. Then the function h^ (J3, k) extended to M. d as a periodic 



function with period 2ir satisfies all requirements of Definition ^. {A with L\ = Li ( [— ir, it 
Proof. The lattice Fourier transform of h (/?, r) equals 

(7.11) 



h(P,k) = e" mk *$(He" imk = ®(P m ) e ~ 

meZ d meZ d 



•im-(k-k») 



Since the above expression naturally defines h(P,k) as a function of k — k*, it is sufficient 
to take k* = 0. To get (|2.34|) . we use the representation of $ (r) in terms of inverse Fourier 
transform ()2.60|) 



$(r) 



(27r) d JRd 



$ (k) e ir - k dk, $(/3m) 



$ [ ^k ] e imk dk. 



(7.12) 
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We split $ (/3m) into two terms: 



* (l3m) = I * (^ { ) 4 G { ) e '"" dC + *' (m) • (713) 

with \l/ (£) defined by (j5.12|) . The first term in (|7.13|) coincides with the inverse lattice Fourier 
transform, its lattice Fourier transform is explicitly given and can be treated as in Lemma 
17.11 The second term gives O (/3 ) with large N for Schwartz functions $. Using these 
observations we check all points of Definition 12.91 as in Lemma 17.11 ■ 



7.3 The Nonlinear Maxwell equation 

We expect that the approximate superposition principle can be generalized to the Nonlinear 
Maxwell equations (NLM) in periodic media studied in j3] . A concise operator form of the 
NLMis 

<9 r U = — -MU + J^nl (U) - J , U (r) = for r < 0, 

Q 

where the excitation current 

J (r) = for t < 0. 

We were studying the properties of nonlinear wave interactions as described by the Nonlinear 
Maxwell equations in series of papers pQ-jEl- Our analysis of the solutions to the NLM uses 
an expansion in terms of orthonormal Floquet-Bloch basis G n ^ (r, k), n — 1, , namely 

oo 

U (k, r, r) = E ^ ( k ' r ) (r, k ) > k e *]* ■ (7-14) 

f=±l n=l 

This expansion is similar to (|2.18|) with J replaced by oo, since the linear Maxwell operator in 
a periodic medium has infinitely many bands. The excitation currents take the form similar 
to forcing term in (|3.1)L namely 

J (r, k, r) = j n , + (k, r) G„„+ (r, k) e "^ (k)T + j n ,_ (k, r) G n ,_ (r, k) e^ T , 

J n (r,k, r) = 0, n 7^ n , 

with a fixed n = uq. The difference with is that time-independent h n ^ (k) is replaced 
by j n £ (k, r). The functions j n ,c (k, r) for every r have the form of wavepackets in the sense 
of Definition 12.91 or in particular the form similar to ()2.36|) with fixed k*. 

The Existence and uniqueness Theorem for the NLM is proven in jl], in particular 
function- analytic representation of the solution as a function of the excitation current. The 
results of this paper can be extended to the NLM equations provided that certain technical 
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difficulties are addressed. Particularly, the classical NLM equation allows for the time disper- 
sion with consequent time-convolution integration in the nonlinear term. This complication 
can be addressed by approximating it with a nonlinearity of the form (|2.22|) with an error 
O (g) = O (/9 2 ), see jH]. Then the derivation of the approximate linear superposition principle 
for wavepackets can be done as in this paper. Another complication with the NLM is that it 
has infinite number of bands. 



7.4 Dissipative terms in the linear part 

Equations ()2.3|) and ()2.6H) involve linear operators iL (k) with purely imaginary spectrum. 
Quite similarly we can consider equations of the form 



<9 T U (k,r) 



-G (k) - -L (k) 

Q 



U(k,r) + F(U) (k,r 



(7.15) 



where a Hermitian matrix G (k) commutes with the Hermitian matrix L (k) and G (k) is non- 
negative. In this case the approximate superposition principle also holds. The proofs are quite 
similar. In the case ([2.61)1 . which corresponds to of PDE, G (k) determines a dissipative term, 
for example G (k) = |k| 2 /,k 6 M. d , where I is the identity matrix, corresponds to Laplace 
operator A. When such a dissipative term is introduced, we can consider nonlinearities 
F which involve derivatives, see jH], [H] in a similar situation. For such nonlinearities our 
framework remains the same, but some statements and proofs have to be modified. We will 
consider this case in a separate paper. 



8 Appendix A: Structure of a composition monomial 
based on oscillatory integrals 

Every composition monomial M \ J-,T, X^X(m)J (^i • • • h m j based on oscillatory integral 

operators J 7 ^ 1 ' as defined by ()3.14|) and the space decomposition as defined by ()5.1|) has the 
following structure. Let T be the tree corresponding to the monomial M. The monomial 
involves integration with respect to time variables T/m where N £ T are the nodes of the 
tree T. The monomial also involves integration with respect to variables k^, N e T. The 
argument of the integral operator M (^F, T, A(§), C(m)) involves only end nodes (of zero rank) 
and has the form 

n h N (k N ) . 

rank(iV)=0 

The kernel of the integral operator involves the composition monomial M (x,T, X(s)X( m )j 
based on the susceptibilities tensors xT? (k, ki m \ \ with the same tree T. Note that the 
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phase matching condition (|3.12j) takes the form 



k N = k' N + . . . + k 



(m(aO) 

N 



a(N)- 



8=1 



Recall that if Cj (iV), i = 1, . . . , /i(iV) is the i-th child node of iV, then the arguments in 
()3.14|) are determined by the formula 

k Ci(Ar) - k w . 



Hence, the kernel of the integral operator M yj 7 , T, X^), C( m )j (hi . . . h m j involves the prod- 
uct of normalized delta functions 

Y[ 5 (k N -k Cl(N) - k CM(JV)( v , j . 

rank(7V)>0 

and the integration with respect to k N is over the torus 



n / )[-] n 



and, obviously, the variable k Nt corresponding to the root node N* is not involved into the 
integration. 

Since every operator at a node iV of the monomial M {? ', T, A(3), C( m )) contains the 
oscillatory factor 



exp < l 



( n uj (k N ) - C n oj (k' N ) 



'(AT) 



A m ), , ( u("0 



T 



(TV) 



we obtain the following total oscillatory factor 



exp < i-$ 



where the phase function $ T ^ ( k ) of the monomial is defined by the formula 



TVeT 



M (/V) 



(>, (k) - £ #<"» w (k Ci{N) ) 



i=l 



T 



(TV)- 



(8.2) 



The vectors k, f and ( are composed of k^, and Ctv using the standard labeling of the 
nodes. 
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Notice then that the oscillatory exponent \8.1\) is the only expression in the composition 
monomial which involves parameter g. Observe also that the FM condition takes here the 
form 

H(N) 
i=l 

The domain of integration with respect to time variables is given in terms of the tree T by 
the following inequalities 

D T = {t (n) : < t (n) < T ip(N)) , NeT\N*} (8.3) 

where p (N) is the parent node of the node N. Using introduced notations we can write the 
action of the monomial M (^F, T, A(|), C(m)) m ^ ne f° rm 

exp j i-$ r>? (k, r) I M (x, T, A, (, k) f[ 

^ ^ ' rank(V)=0 

JJ 5 (k N - k Cl(JV) k CKN)(Ar) J JJ dk w JJ dr (JV) . 

rank(2V)>0 N^N* N^N, 

Note that m equals the number of end nodes, that is nodes with zero rank and they are 
numerated using the standard labeling of the nodes, that is 

hi (ki) • • • h m (km) = II hjv ( kiV ) • 

rank(7V)=0 

The formula (|8.4|) gives a closed form of a composition monomial based on oscillatory integral 
operators with an arbitrary large rank. 



9 Appendix B: Proof of the refined implicit function 
theorem 



Here we give the proof of Theorem 14.251 

First, we consider the following elementary problem which provides majorants for the 
problem of interest. Let a function of one complex variable be defined by the formula 



f{u)=C T Y, umR r m = CT 



m=2 



u 2 /R\ 
1 - u/R r _ 



T 



Cf>0,Rr> 0. 



(9.1) 
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In this case {x\ . . . x m ) = CfRj^xi . . . x m . Let us introduce the equation 

u = T (u) + x, u, x G C (9.2) 

which is a particular case of (J4.1j) . A small solution u (x) of this equation such that u (0) = 
is given by the series 

oo 

u = G (x) = J2 G {m) x m , 

m=l 

which is a particular case of formula (|4.14|) . The terms G^x m of this problem are determined 
from (|4.18jl and can be written in the form (|4.29j) 

G (m) x m = c T M (JF, T) x m . (9.3) 

Obviously, 

M (JF, f) x m = C i PR~ e{T) x m (9.4) 

where i (T) is the incidence number of the tree T, e (T) is the number of edges of T. Now 
we compare solution of the general equation (|4.1j) . It is given by the formula (J4.14j) with 
operators (u m ) admitting expansion (J4.29|) . Since 

||jrM|| < C T R7 T m , 
where the constants are the same as in (J9.1)) we have 

\\M (F, T) ( Xl . . . x,) || < M (F, T) || Xl || . . . ||x,|| , 

implying 

c T ||M(^,T)( Xl ...x m )|| < ctM(F,T) ||x 1 ||...||x m || =G( m )||x 1 ||...||x m ||. 



TeT m TGT„ 



Solving (|9.2|) we get explicitly 



(9.5) 



-=t( i -v rr ^5) =e(i) ' ,!= S +i - 



We have the following estimate of the coefficients 



G {m) < WTT^-B ^ ( " , m = 1, 2, . . . , (9.6) 



2{C T + R r )y R% , 
(see |1] for details in a similar situation). From ()9.4j) and ()9.6|) we infer the following inequality 



rei) 



which hods for all CV, R T > 0. We set = R T = 1 and obtain the desired bound (|4.35j) . 
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10 Notations and abbreviations 



For reader's convenience we provide below a list of notations and abbreviations used in this 
paper. 

AFM alternatively frequency matched, see Definition 15. 191 

ANFM - alternatively non- frequency- matched see Definition 15.191 

band-crossing points see Definition 12 .HI 

cc - complex conjugate to the preceding terms in the formula 

composition monomial — see Definition 14.91 

decoration projections see f!4.36|) . (|4.37JI 

decorated monomial see Definition I4.2UI 

CI monomials - cross-interacting monomials, see Definition 15.81 

FPU, Fermi-Pasta-Ulam equation - (jUUD , (l2~TTl) . fl7~TJ) 

Floquet-Bloch modal decomposition — see (I7.14J1 

Fourier transform see (I2.59J) 

FM - frequency matched, see Definition I5.1UI see also (J5.42|) 
homogeneity index of a monomial Definition 14.91 
homogeneity index of a tree Definition 14. Ill 

incidence number of a monomial — number of occurrences of operators in the 
composition monomial 

incidence number of a monomial — see Definition 14. 1UI 
incidence number of a tree — Definition 14.121 
lattice Fourier transform see f!2.2|) 
monomial - Definition 14.91 

NFM - non-frequency-matched see Definition 15. 101 see also ()5.46|) 
oscillatory integral operator see (|3.8|) . (|3.3|) 
rank of monomial see Definition 14.91 
root operator (|4.20j) 

SI monomials - self-interacting monomials, see Definition 15.81 

Schwartz functions — infinitely smooth functions on IR d which decay faster than any 
power, see ()7.6|) 

single- mode wavepacket — see Definition 12.91 
submonomial (|4.1U|) 
wavepacket see Definition | 



dim ~ 1)k = (^)CT dk ' • • • dk(m_1) - see ra> 

(m-l)d _ Miy am p m(m-l)d 



-7T, 7r - see 



$FI^ or D m = see (gH^ 



E = C([0,t*],L 1 ) see (j23H|) 

p{m) _ m _ij near p era tor in L 1; see (|2~2ll . (l2~T)4l 

■p( m ) _ Das i s element of the m-linear operator in E see (13.81) 

F^l see Km 
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/i£ (/?, k) , ( = ± - Fourier transform of the wavepacket initial data (j3, r), see Definition 
( s^J ' C = =t ~ Fourier transform of the wavepacket (/3r)initial data, see Definition 

h* (k, /3) - a function nullified outside /3 1_e vicinity of ±k*, see (|5.15j) 

k = (ki, . . . , kd) G [— 7r, 7r] - quasimomentum (wave vector) variable, (|2.2j) . (|2.25|) . 

k = (fei, . . . , k d ) e M. d - Fourier wave vector variable, (|235)) . (|2~25|) . 

k* = (A;*i, . . . , - center of the wavepacket see Definition 12.91 

\t*i - center of /-th wavepacket 

k = (k', . . . , k^ m ^), - interaction multi-wave vector, (|2.25J1 . ()3.7|) . 
kW (k,k) =k-k' - k^" 1 ) - see ()2~23D 

L 1 - Lebesgue space L x {\— it, 7r] d J or L x (R d ) - see ()2.31)) and ()2.66j) 
n - band number 

n = (n', . . . , n^) - band interaction index, ()3.7|) 
Vr = " ' > £2) ~ s P atial gradient 

O (fi) - any quantity having the property that is bounded as /i — > 0. 

(k) = £u; n (k) - dispersion relation of the band (C,^), see (12.13)) 
u' nQ (k) = V k ^ no (k) - group velocity vecto r 

u n (k)- n-th eigenvalue of L (k), see (J2.13)) : dispersion relation of n-th band 
^ - cutoff function in quasimomentum domain, see 1)5.12)1 

(f), H {k, kj = (u n (k) —('u n i (k') — . . . — (^u n ( m ) (k^) - interaction phase function, ()3.9j) 
7r - see (15.13)) 

II n> £ (k) - projection in C 2J onto direction of g n ^ (k); see (J2.19)) 
r = (n, . . . , Td) - spatial variable 
g = (3 2 - (J23I3) 

a - the set of band-crossing points, see Definition 12.31 
U (k) - Fourier transform of U (r), see (12.59)) 

U n>c (k, r) = u U)C (k, r) e _i F Cwn(k) - amplitudes, see (f3~2|) 
£ = ±or£ = ±l- band binary index. 

£ = k', • • • , ~~ binary band index vector, see ()3.7)1 
Z* - complex conjugate to Z 
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